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SOLUTION OF AN INITIAL-VALUE PROBLEM 


IN LINEAR TRANSPORT THEORY: MONOENERGETIC NEUTRONS 
IN A SLAB WITH INFINITE REFLECTORS* 

By Perry A. Newman and Robert L. Bowden** 

Langley Research Center 

SUMMARY 

The solution of an initial -value problem in linear transport theory is obtained by 
using the normal-mode expansion technique of Case. The problem is that of mono- 
energetic neutrons migrating in a thin slab surrounded by infinitely thick reflectors and 
the scattering is taken to be isotropic. The results obtained indicate that the reflector 
may give rise to a branch-cut integral term typical of a semi- infinite medium whereas 
the central slab may contribute a summation over discrete residue terms. Exact 
expressions are obtained for these discrete time eigenvalues, and numerical results 
showing the behavior of real time eigenvalues as a function of the material properties of 
the slab and reflector are presented. These eigenvalues are finite in number and may 
disappear into the branch cut or continuum as the material properties are varied; such 
disappearing eigenvalues correspond to exponentially time-decaying modes. The two 
largest eigenvalues can be compared with critical dimensions of slabs and spheres, and 
the numerical values are shown to agree with the criticality results of others. In the 
limit of purely absorbing reflectors or a bare slab, the present solution has the same 
properties as have been previously reported by others who used the approach of Lehner 
and Wing. 


INTRODUCTION 

Linearized transport equations are encountered in a number of different areas such 
as neutron diffusion, radiative transfer, sound propagation, and plasma theory, and the 
extent to which they correspond to reality varies from problem to problem even within a 


The material presented herein was a thesis entitled "Time Dependent Mono- 
energetic Neutron Transport in a Finite Slab With Infinite Reflectors" submitted in partial 
fulfillment of the requirements for the degree of Doctor of Philosophy in Physics, Virginia 
Polytechnic Institute, Blacksburg, Virginia, December 1969 by Perry A. Newman. 

**Associate Professor, Virginia Polytechnic Institute and State University. 




given area. Since there are some problems where the approximations required for 
linearization are not unduly severe from a physical point of view, a substantial effort to 
develop methods for solving such equations has been made. However, analytical solu- 
tions have been obtained for only a small class of highly idealized problems and few of 
these solutions have been for the time -dependent situation. Therefore, it might be 
expected that any new time -dependent result should provide additional insight into the 
general character of such solutions. In any case, analytical solutions provide a set of 
check cases for comparison of the approximate methods (see, for example, refs. 1 and 2) 
which are used in practical applications. The present idealized neutron transport prob- 
lem, being somewhat more complicated than those considered previously, gives such new 
results. In addition, the procedure used to obtain the present solution is generally appli- 
cable to other time-dependent problems in linear transport theory and is therefore of 
interest in its own right. 

The basic assumptions customarily made in neutron transport theory (see, for 
example, refs. 1, 2, and 3) have been summarized by Wigner (ref. 4). Since the conven- 
tional derivations of the governing transport equation are based mainly on plausibility 
arguments, Osborn and Yip (ref. 5) have reexamined the situation starting from a micro- 
scopic point of view using quantum mechanics.* They conclude that they are only par- 
tially successful in justifying the conventional neutron transport equation because a num- 
ber of required approximations are merely stated and not analyzed or evaluated; however, 
their approach "brings many aspects of the neutron problem into contact with other micro- 
scopic transport theories." For the present problem, the additional restrictions to mono- 
energetic (called one-group or constant cross-section approximation in ref. 1) neutrons, 
plane geometry, and isotropic scattering in the laboratory system are made. (See refs. 1 
to 4.) It is pointed out in references 1 and 2 that even though the monoenergetic approxi- 
mation is rather severe from a physical point of view and is made primarily in order to 
obtain analytical solutions which usually cannot be obtained for the general case, it forms 
the basis of the more physical multigroup approximations. Under these restrictions, 
the solution of the initial-value monoenergetic neutron-transport equation for a one- 
dimensional slab of finite thickness surrounded by infinitely thick reflectors of a different 
material is obtained in this report. 

One mathematically rigorous approach which has been used to treat such problems 
is a spectral analysis. Lehner and Wing (refs. 6 and 7) used this approach to solve the 
initial value, monoenergetic neutron-transport problem for a bare slab where the 

*Osborn and Yip (ref. 5) give four reasons for using a quantum mechanical treat- 
ment: (1) a formalism for describing the creation and destruction of particles exists, 

(2) the neutron-nuclear interactions are truly quantum phenomena, (3) the interpretation 
of an observable density in phase space, and (4) other peculiarly quantum effects such as 
spin and associated statistics. 


2 



scattering is isotropic. Lehner (ref. 8) considered a slab of finite thickness surrounded 
by a purely absorbing medium which had the same total macroscopic cross section as the 
slab. Very recently Hintz (ref. 9) has generalized this problem by allowing the purely 
absorbing medium to have any cross section. Mika (ref. 10) has studied the initial -value 
problem for monoenergetic neutrons in a nonuniform slab surrounded by a vacuum but did 
not obtain results as complete as those of Lehner and Wing (refs. 6 and 7). In particular, 
theorems concerning the reality and number of discrete time eigenvalues were not 
established. 

Another approach which has been used to solve a few time -dependent, monoenergetic 
neutron-transport problems in plane geometry is the normal-mode expansion technique of 
Case. (See refs. 2 and 11.) This method was used by Bowden and Williams (refs. 12 
and 13) to analyze the problem which had been treated by Lehner and Wing (refs. 6 and 7). 
A second application of this technique was made by KugCer and Zweifel (ref. 14) to the 
time -dependent, one-speed albedo problem for a semi-infinite medium. Finally, Erdmann 
and Lurie (refs. 15 and 16) have also utilized this approach in a two-media time -dependent 
problem, the time decay of a plane isotropic burst of monoenergetic neutrons introduced 
at the interface of two dissimilar semi-infinite media. In all these time-dependent solu- 
tions, contributions due to various parts of the spectrum of the transport operator have 
been indicated by suitably deforming the integration contour of the inverse Laplace 
transformation. In view of these successful applications of Case's technique, in particu- 
lar references 15 and 16, this technique has been chosen to analyze the present problem. 

In this problem, one would expect to find discrete time eigenvalues (time constants) and 
obtain some insight concerning their behavior as a function of material properties. Since 
the reflectors can scatter as well as absorb neutrons, the solutions for the bare slab and 
slab surrounded by purely absorbing media are included and it is shown that the present 
solution agrees with such solutions (refs. 6 to 9) for these special cases. Some pre- 
liminary results of the present work were given in reference 17 and a summary of the 
present results is given in reference 18. 


SYMBOLS 


A 


nondimensional slab half -thickness (see eq. (118)) 


A m , a m>fr m expansion coefficients in Case's normal-mode expansion for medium m 
(see eq. (29)) 

A m± , a m ± definite parity expansion coefficients (see eq. (32)) 
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®m±»^l± 


C* 

c 

c 


m 


m± 


'm 


E 


m± 


F m±> F ± 


f m± 

^m 

h m 

Im( ),Re( ) 


m± 


’m± 


m± 


slab half -thickness 

expansion coefficients for solution of associated eigenvalue problem 
(see eq. (76)) 

contour in z’ -plane (see fig. 3) 
contours in s-plane (see figs. 2 and 4) 

expansion coefficients in full-range normal-mode expansion of initial 
distribution in medium m (see eq. (47)) 

contour in s-plane (see fig. 4) 

mean number of secondary neutrons per collision in medium m 

expansion coefficients defined by equations (50) 

integrations over initial distribution f(x,/i) given by equation (46) 

initial neutron angular flux, generally referred to as the initial 
distribution 

definite parity parts of initial distribution in medium m (see eq. (16)) 
solution of equation (42) 
given by equation (80) 
imaginary and real parts 

inhomogeneous terms given by equations (55) and (57) 
inhomogeneous terms given by equations (56) and (58) 
given by equation (62) 

integrations over initial distribution f (x,/i) defined by equation (67) 
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lower limits given by equation (68) 


'm 


M, 


integrations over initial distribution f (x,p.) given by equations (74) 
and (75) 


A ’m 

P 

^me>^mi 

s 

s n 

t 

^0>^0m>^m 

x 

z 

a l±> a 2± 

h±>P2± 

y 

«m( s ) 

5(v - ji) 


given by equation (80) 

denotes that Cauchy principal value is to be taken upon integration 

denote regions in s -plane (see figs. 2 and 4) 

complex variable of Laplace transformation (see eq. (6)) 

denote values of s for which associated eigenvalue problem has 
nontrivial solutions 

real time multiplied by constant neutron speed 

Case's X-functions given by equations (59) to (61), respectively 

geometric coordinate perpendicular to slab faces (see fig. 1) 

complex variable 

defined in equation (28) 

given by equation (89) 

given by equation (90) 

associated with inverse Laplace transformation (see eq. (7)) 
defined by equation (30) 

Dirac delta function 




given by equation (118) 



£ n denote values of £ for which associated eigenvalue problem has non- 

trivial solutions 

£g £ n with largest real part 

A. m defined by equation (25) 

p. direction cosine (see fig. 1) 

v introduced as complex separation parameter in equation (22) but used 

thereafter as real, z being used to denote complex values 

±^Q m values of z for which O m (z,s) = 0 

a m total macroscopic cross section for medium m 

°R> a D given by equation (118) 

a min = min (°i>°2) 

<$ given by equation (110) 

<P mv continuum mode in Case's method (see eq. (24)) 

cp discrete mode in Case’s method (see eq. (26)) 

±!/ 0m 

41 neutron angular flux 

i {/ definite parity parts of Laplace transform of neutron angular flux, 

generally referred to as transformed solution 

4 / m± 4 / ± i n medium m (see eq. (17)) 

Vc±) ’/'mpi P arts of ’/'mi ( see e( l* ( 31 )) 

an elementary solution (see eq. (22)) 

solution of associated eigenvalue problem 
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S2 m dispersion function (given by eq. (27)) 

12^ given by equation (49) 

Subscripts: 

c complimentary solution 

e,i regions exterior and interior to a curve 

m physical medium, 1 for reflector and 2 for slab 

n eigenvalue, 0 for the one with largest real part 

p particular solution 

u solution due to part of initial distribution which has not been scattered, 

that is, uncollided 

v associated with continuum modes 

± definite parity, + for even and - for odd 

Superscript: 

± limiting values of a function on its branch cut as argument approaches 

cut from upper (+) and lower (-) half-planes. 

A bar over a symbol denotes the associated eigenvalue problem quantity. 

PROBLEM DEFINITION 

Basic Equations 

Consider a slab of material which scatters neutrons isotropically (in the laboratory 
system), extends from x = -a to x = a, and is characterized by a total macroscopic 
cross section 03 and a mean number of secondary neutrons per collision C2. This 
uniform slab is surrounded by uniform infinitely thick reflectors of another material 
characterized by the nuclear properties and c^. (See fig. 1.) For isotropic 
scattering of monoenergetic neutrons in a sourceless medium with plane geometry, the 


7 




Reflector Slab Reflector 

Figure 1.- Geometry of problem. 


neutron angular flux ^(x,/u } t) satisfies the equation (see refs. 1, 2, and 3 for a complete 
list of the approximations and assumptions) : 


a*(x,M,t) , „ 8*(x,/i,t) 


at 


+ JLt- 


9x 


+ a(x) *(x,/i,t) = 


c(x) ct(x) 


f 1 *(x,p',t) 
J_1 


d/i' 


(1) 


where t is the real time multiplied by the constant neutron speed, x and p. are 
shown in figure 1, and a(x) and c(x) are given by 


^1> C 1 

o(x),c(x) =/ 

(|x| > a. 

L°2’ c 2 

(|x| <a 


Since multiplying media are not of infinite extent, Cj has been taken to be less than 
unity. Therefore, equation (1) is to be solved subject to the boundary conditions 

lim *(x,p,t) = 0 ( | x | -*• (3) 

and the continuity conditions 

*(±a+,M,t) = ^(±a-,/i,t) (4) 
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given the initial condition 


*(x,/i,0) = f(x,M) (5) 

which is assumed to satisfy equation (3) and be extendable without poles or branch cuts in 
the finite p. -plane except perhaps for a discontinuity across the imaginary axis. When 
the material properties of the reflectors are taken to be those of a vacuum, this problem 
reduces to that of references 6 and 7 whereas for purely absorbing medium, it reduces to 
that considered in references 8 and 9. When a(x) is constant and c(x) = 1 (that is, 
neutrons are conserved) equation (1) is the one -dimensional Boltzmann equation for the 
Lorentz model of kinetic theory (ref. 19). 

The general procedure used here to solve the mathematical problem presented in 
equations (1) to (5) consists of the following steps: 

(1) Remove the t-dependence with a Laplace transformation 

(2) Solve the transformed equation by applying Case's technique 

(3) Determine the analytic properties of this transformed solution in some right- 

half s -plane 

(4) Recover the t-dependence and simplify by suitably deforming the integration 

path of the inverse transformation 

(5) Calculate real discrete time eigenvalues as a function of material properties if 

and when they exist. 

Since many details are involved in performing these few steps, a brief synopsis is 
given. Step (1) is easy and the transformed equation and boundary conditions are given by 
equations (8) to (10). Step (2) is accomplished by construction of a solution from Case's 
elementary solutions, or normal modes. Since the geometrical symmetry aids in this 
construction, symmetry properties are introduced immediately after the time removal. 
The elementary solutions of Case are given by equations (24) to (28). For details con- 
cerning these solutions and their completeness and orthogonality properties, the reader 
is referred to references 2, 11 to 16, and 20. Appendixes A and C summarize the 
important results taken from these references which are required in this report. Con- 
struction of the transformed solution is done in the section so titled and a few details are 
given in appendixes B and D. Equations (50) to (58) give the expansion coefficients of the 
transformed solution implicitly. Extension of these equations to the complex plane is 
presented in appendix E. Since the solution valid for all regions of the transform plane 
is needed for step (3), it is obtained in the same section and is given by equations (70) to 
(73). Step (3) is performed in the section "Properties of Transformed Solution" and 
many of the details are given in appendixes F, G, and H. In particular, one must examine 
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where the associated eigenvalue problem has nontrivial solutions and show how these 
solutions enter the transformed solution. Details for step (4) are found in the section 
"Recovery of Time -Dependent Solution" and appendix I. The previously cited results 
(refs. 12 to 16) lead one to expect that the reflectors should contribute continuous- 
spectrum type terms typical of a semi-infinite medium whereas the central slab should 
give rise to some point-spectrum type terms and their corresponding discrete time 
eigenvalues. The solution ^(x^t) is given by equation (114). Step (5), the calculation 
of real time eigenvalues, is outlined rather explicitly in appendix J and the numerical 
results are presented and discussed in the section "Calculation of Time Eigenvalues." 

The report is concluded with a short section showing how for special values of the nuclear 
properties the present solution reduces to those obtained previously by others (refs. 6 
to 9) who used a different method. 


Time Removal 

If one takes the Laplace transformation of 'Ji r (x,ji,t) as 

^(x,m,s) = f e" st *(x,ju,t) dt (6) 

J 0 

then the inverse transformation required to recover the t-dependence is 

1 py+ico 

tf(x,/i,t)=— \ e st !//(x,ju,s) ds (7) 

2lU dy-ioo 

where y is to the right of all singularities and branch cuts of i//(x,/i,s) in the trans- 
form plane (s-plane). From previously cited work of others, it is expected that the path 
of integration in equation (7) can be deformed to indicate more precisely the character 
of ^(Xj/ijt). When the transformation of equation (6) is applied to equation (1), integra- 
tion by parts is performed in the usual manner, and use of the initial condition (5) is 
made, the following expression is obtained: 


ju^^(x,jli,s) + (s + a(x)]^(x,^,s) = — X v,- — j* ^(x,M’,s) d p.' + f (x,m) (8) 


Equations (3) and (4) become under the same transformation 


lim i//(x,|Lt,s) = 0 


(|x| - °°) (9) 
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and 


^(±a+,M,s) = if/(± a-,Ai,s) (10) 

Symmetry Considerations 

Before applying Case's technique to solve equation (8) subject to conditions (9) 
and (10), it is useful to examine some symmetry properties of the transformed solution 
which follow directly from the governing equations. In references 12 and 13, these ideas 
were introduced at a later step, but here they aid in the construction of the solution. An 
arbitrary function of two variables g ( x,jjl) can be written as the sum of its even and odd 
parts, namely, g + (x,p.) and g_(x,/i.). They are given, of course, by 

g ± (x,/x) =|[g(x,M) ± g(-x,-pj] (11) 

and have the property 


g ± (-x,-jn) = ±g ± (x,ju) (12) 

Since c(x) and ct(x) are even functions of x, it is easily shown from equation (8) that 
the even and odd parts of i//(x,^,s) obey the equation 

M-^^ ± (x,M,s) + [s + a(x)]^ ± (x,ji,s) = J i dn' + f ± (x,/i) (13) 

The boundary conditions for i// ± corresponding to equations (9) and (10) are written as 

lim \p ± (x,ix,s) = 0 ( I x | - «°) (14) 


and 


^ ± (a+,M,s) = i//±(a-,/i,s) (15) 

where the ± subscripts denote definite parity parts of a function, that is, even and odd. 
(See eqs. (11) and (12).) Equations (13) to (15) indicate the following: 

(1) All solutions of the homogeneous equation associated with equation (13) can be 
made to have a definite parity. 

(2) The boundary conditions preserve the parity. 
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(3) The definite parity parts of an initial distribution excite inhomogeneous solu- 
tions of corresponding definite parity. 

Therefore, this problem can be separated into two problems, one for i// , the other 
for j p_, and the results can be combined at any stage of the calculation. The functions 
f ± (x,p) and i// ± (x,p,s) are broken up as 


f ± (x,p) 


fl±(x,p) 

f2±( x ,M) 


and 


rf/ ± (x,n,s) 


s) 

i// 2± (x,p,s) 


(|x| > a) 
(|x| < a) 


(16) 


(|x| > a) 
(|x| < a) 


(17) 


so that equations (13) to (15) become 
9 


cr w 
m m 


y ^mi^M'jS) dp' + f m± (x,p) (18) 


^^m±( x ^’ s > + ( s + a m)^ m ±( x ’^ s ) = 2 

where throughout the subscript m = 1 denotes medium 1 and m = 2 denotes medium 2, 


lim i// 1± (x,p,s) = 0 


(|x| - °°) 


(19) 


and 


i// 1± (a,p,s) = ^ 2± (a,p,s) (20) 

The notation g m± (a,p) means the limit of g ± (x,p) as x - a from medium m. 

Elementary Solutions 

Solutions of equations (18) are constructed from Case's elementary solutions which 
are denoted here as i// mj ,(x,p,s). These elementary solutions are solutions of the homo- 
geneous equation corresponding to equation (18); that is, 

+ ( s + °m)tfWk»^ s > = 2 c m a m J ,*) d#i' (21) 

in the form 

= <Pm iXm-,®) e“( S+a m) x /^ (22) 
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where v is a complex parameter introduced in this separation of variables and 
< Pm l /(J Lt J S ) i- s normalized as 

f <P m v (^,s) dp = s + o m (23) 

J -1 

These solutions have been investigated in references 12 to 16 and many of their results 
are given in appendix A and are used herein. They show that the values of v for which 
solutions <p m j,(p,s) can be found are v real (-1 = f = 1) and for some region of the 
s-plane v = ±Pg m . For -1 = ^ = 1, the solutions are 

cp mv {li,s) =|c m a m ^P 17 ^+ X m {v,s) - p) (24) 

where P denotes that the Cauchy principal value is to be taken upon integration, 

6(v - p) is the Dirac delta function, and X m (^,s) is determined from the normalization 
as 


VK 3 ) = s + u m - c m a m v tanh' 1 v 


(25) 


These are called the continuum modes and exist for all values of s. There are two 
discrete solutions 


1 c m a m z '0m 

<P ± i (M,s)=^— — 

Om 2 t'Oni + P 

at v = ±^Q m provided that the function S2 m (z,s) 


^m( z J s ) = s + ^m “ c m°m z tanh 


■li 


(26) 


(27) 


of two complex variables s and z vanishes at the two points z = ±t'Q m (s). This con- 
dition occurs when s lies inside the curve C m (s e S m ^, see fig. 2) defined by (see 
refs. 12 and 13) 


Cwi = < 


s + °m 
'm u m 


= a' + i/ 3* 


a .^JL Unh -i( 2 JL 


(28) 


Note that i^Q m is an analytic function of s for s e S m ^ except for a branch cut on the 
real s -axis between -CT m and _CT m (l - c m ) and that +^om denotes that zero of 
fi m (z,s) for which > 0 when Re(s) > - a m (l - c m ). The important result is 

that the general solution of equation (21) can be expressed as the linear combination 
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iirc a / 2 
m m 


Re(s) 


-irrc a 
in 


J 2 


Figure 2.- Regions in a single-medium s-plane. Location of Im(s) axis 

depends on whether ^ 1. 


where 5 m (s) is defined as 
6 m (s) h 




1 

10 


5 m (s) + f A m (f) ^ m „(x,p,s) dv 
J -1 


( s e S mi) 
( s e s me) 


(29) 


(30) 


and the s-dependence of the expansion coefficients has not been determined. Note that 
the present notation is slightly different than that used by other authors. 

CONSTRUCTION OF TRANSFORMED SOLUTION 

Solutions of equation (18) are now obtained by constructing even and odd particular 
solutions i// m p ± (x,;u.,s) and adding to them solutions of the corresponding homogeneous 
equations i^ mc± (x,q.,s) so that conditions (19) and (20) can be satisfied; that is, 

xl/ m± {x,H,s) = ^ mc± (x,M,s) + ^ mp± (x,p,s) (m = 1, 2) (31) 

14 



The functions 'P mc± and ^ m p ± are constructed from Case's elementary solutions 

One must select the expansion coefficients in a general expansion, such as 
equation (29), so that the given boundary conditions are satisfied. 


Explicit Form of & mc± 

The solution in the form of equation (29) does not have definite parity. However, 
for a medium which is connected and symmetric about x = 0 (such as the slab), even 
and odd solutions ^2c± can be written as 


^2c±( x ^> s ) 


= a 2±fyv n J x >^> s ) ± (x,M,sj|5 2 (s) 

_ U Li (J Cl _ 


y A 2± ( v) j~^ 2 v (x, M , s) ± ^2 ( _ v ) (x, M , s jj dv 


(32) 


where the expansion coefficients have been redefined as 

a 2 ± = |( a 2 ± b 2 ) 

A 2±( y ) =|[A2(^) ± A 2 (-i>)j 

J 

Note that the properties 

^W" X '~^ s)= ^0m (X ’^ S)> 


(33) 


(34) 


have been used. 

For a medium which extends to infinity in the x-direction (such as the reflectors), 
the boundary conditions (9) require that 

b x = A 1 (-i/) = 0 (0^1/51 if x-H (35a) 

or 

a t = A 1 (y) s 0 (0=^ = 1 if x — -°°) (35b) 
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for the expansion coefficients in equation (29) when Re(s) > -crj. Thus, in the reflectors, 
the solution has the form 


^ 1 ( x ,/ j ., s ) =/ 


bi i//_^(x,p,s) 5 1 (s) + J A t (y) ils lv (x,ii,s) dv (x < -a) 

pi y 

ai > 2 , (x,fi, s ) 5 1 (s) + ^ A 1 , (V) i^ lz ,(x,p,s) dv (x > a) 


(36) 


for Re(s) > -a-j. . The continuity conditions (10) and the parity of the solutions ^2c± are 
used next to relate the primed and unprimed coefficients in equation (36). One finds that 
an even solution inside the slab requires 


ai' = bi 
Ai(p) = A^-v) 


(0 S v $ 1) 


(37a) 


whereas an odd solution inside the slab requires 

a i' = -bl 
A{(v) = -A 1 (-i') 

In view of equations (37a) and (37b), one defines, respectively, 


a l+ = b l 


A 1+ (-^) = Ai(-^) 


and 


a l- £ 

A 1 _(-v) = A^-v) 



(0 £ v s l) 


(37b) 



(0 ^ s 1) 


(38a) 


(0 ^ v S il 


(38b) 


In terms of these coefficients, ^ c± can be written from equation (36) as 
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lp lc± (x,ll,s) = l 


a l±^-^ 01 ^ x ’ M ’ s ) 5 l( s ) + J Q A i±(" y ) ^(.^(Xj/AjS) dv (x < -a) 

±a l ±^u n M’^ s ) M s ) ± f A i ± (“^) rl/ lv k,H,a) dp (x > a) 
ui 


(39) 


for Re(s) > -a^. 


Explicit Form of i^ m p ± 

The definite parity particular solutions are constructed from two particular 

solutions !^ m p as 




mp± 



(x,M, s) ± ^ mp 


( 


-x,-n,s 


'] 


(40) 


The solutions ftmp are obtained in a conventional way by integration of a Green's func- 
tion g m for medium m over all the medium; that is, as 




mp( x ’^-> s ) ~ C 


Medium m 


g m (x,M;x 0 ) dx Q 


(41) 


The function g m satisfies the equation: 



g m (x,M;x 0 ) + (s + o m )g m (x,M;x 0 ) = 


c m a m P 1 
2 J_! 


g m ( x ^ , ; x o) d ^' 


+ 6 ( x - *0) f m( x 0^) 


(42) 


This equation is seen to be the homogeneous equation corresponding to equation (18) for 
x/xg. When one integrates over all xq in medium m, it is seen that equation (18) is 
obtained. Note that g m is not exactly what is customarily called the Green's function, 
since 6(x - xq) has been weighted with f m (xg,p.). Upon integrating equation (42) on x 
from Xq - e to x 0 + e and taking the limit e — 0, one obtains the jump condition 

gm( x 0+>M;x 0 ) - g m (x 0 -,/r,XQ) = (43) 

In appendix B, g m is constructed from Case's elementary solutions and, as a result, the 
explicit forms of ^2p± and ^lp± can de written as 
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^2p ±( x >^> s ) = |^2±( x ^02’ s ) ^i/ 02 (x ’ m ’ s) * F 2 ±( -x ’ i/ 02’ s ) ^-j/ 02 ( x ^,sjj6 2 ( s ) 

+ J 0 f 2±( x >^ s ) ^2 y( x >M,s) di/ 

± F 2± (-x,*,s) ^ 2 (_j,)(x,fi,s) di^ 


^lp±( x ’^’ s ^ = 


F l±( x ^01> s ) ^ 01 ( X .M, s) 

+ [Fi ± (x,-^ 0 1 ,s) - F ± ^a,^oi,s) >/'- I , 01 (x,M,s)J S^s) 

+ f F 1± (x,z/,s) xf/ lp (x,ii, s) d^ 

J 0 

+ jF 1± (x,-v,s) - F ± (-a,i/,s)j ^i(_j,)(x,m,s) dt- (x < -a) (45a) 


for Re(s) > -cr^, and 


^l p± (x,/i,s) = |±[-F ± (-a,v 01 ,s) + F 1 ± (-x,-i/ 01 ,s)]^ o 1 (x,m,s) 


± F^^x^oi.s) 4/_ v (x,ij.,s)\ 5 1 (s) 


± ^ ^F ± (-a,^,s) + F 1± (-x,-i',s)]^ ll/ (x,M,s) 
± J Q F 1± (-x,^,s) ^i(_ y )(x,M,s) d^ 


(x > a) (45b) 


for Re(s) > -a^, where 


F ± (-a,o>,s) = F 1± (-a,-a),s) + F 1± (-a,o>,s) 

C x t \ (s+0'o')x n /o> 

F 2± (x,co,s) = j C 2±( x 0> w ) eV ^ dXQ 


r x / \ (s+a 1 ')x n /w 

F 1± (x,cj,s) = j c l±( x 0> w ) e 7 dx c 
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where u> = and v, -1 £ v ~ 1. Here the C m± are full-range expansion coefficients* 
of the function f m± (x,p,)/p. and are given by 


c m±( x 0^) 


1 

^n>,s) S2 m "(j/,s) 



f m±( x 0’^) 


and if s e S 


mu 


^m±( x 0^0m) “ jn ' ( \ f i dp. 

c m%^0m 2 fi mhm> s ) J “ 1 0m 

^m±( x 0’ -y 0m) = 2771 n J 1 f m±( x 0^) <P- V ,(^) d M 

c m % v 0m a m (- I '0ra> s ) J -1 0m 


) (47) 


Throughout, the + and - superscripts are used to denote the limiting values of a func- 
tion on its branch cut as the argument approaches the cut from the upper (+) and lower (-) 
half-planes. The function £2 m (z,s) of equation (27) has a branch cut along the real 
z-axis (-1,1) where its limiting values are given by 

(^s) = ^ m K s ) ± (-1 £ 1) (48) 

The functions ^^(ZjS) are defined by 

^m( z > s ) = “j~^m( z > s ) (49) 


Equations for Expansion Coefficients 

Solutions in medium 1, |x| > a, have been constructed so that the boundary condi- 

tion (19) is satisfied. Application of the continuity condition (20) permits the determina- 
tion of the unknown expansion coefficients of 4 / mc± which are implicit in equation (31); 
that is, if one substitutes x = a in equation (31), applies the continuity condition (20), 
and uses the explicit forms of i// mc± given by equations (32) and (39), a two-media full- 
range expansion involving the (p mv which contains unknown coefficients a m± and A m± 


*Note that the parity of these coefficients is opposite that indicated by the ± sub- 
script because of the 1/p factor. Nevertheless, the solutions ^ m n± are easily seen to 
have the indicated parity. y 
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is obtained. The same expansion is, of course, obtained for x = -a. This type of 
expansion and its orthogonality relations are discussed in appendix C and it is shown in 
appendix D that such an expansion is obtained for the present problem. Erdmann (ref. 15) 
proved completeness theorems which apply in such time -dependent problems while 
Kuscer, McCormick, and Summerfield (ref. 20) derived orthogonality relations which are 
applicable to two-media expansions which arise in time -independent problems. In 
appendix C, their results are extended to obtain orthogonality relations which are valid 
for all regions of the transform plane. As usual in problems involving a slab, one cannot 
obtain closed-form solutions for the expansion coefficients. However, the orthogonality 
relations (appendix C) can be used to obtain expressions which give the expansion coeffi- 
cients implicitly; that is, the continuum coefficients A 2± (v) are given as the solutions of 
Fredholm integral equations and all the other coefficients are obtained from A 2± (i/). For 
later convenience, however, expressions are obtained for these coefficients in the form 


E2±(2',s) s A 2± (v) ^2"(^,s) €1 2 {v,s) e( S+<72 ) a /^ 
El±( y ,s) s A 1± (-i/) J2j>,s) n{(v,s) e ( s+a l) a A 


(50) 


Use of the orthogonality relations leads, after some algebra, to the following equations: 


E 2± (v,s) = I 2± (v) ± 


k fi 2(°°> s ) v / x 


r i 

Jo E2> ’ s) -^ 


-2(s«2)a/M Xo( .^ s)lx d(1 


2 (M,s) S2 2 (m, s) (m + v) 


6 2 (s) a 2± e "( s+CT 2) a / z '02 x Q (-v 02 ,s) u~T 


(0 £ v ^ 1) (51) 


| c 2 a 2^02^2( l; 02» s ) a 2± e ^ 2 ) a / 02 


-r / \ fk ^2(°°> s ) v / , e 2 ( s+g 2) a /^ X 0 (-m,s) m d # 

2±(^02) ± [2 O^s) °( I ' 02 ’ S )jj 0 2± M,S S 2 2 (m,s) (m + y 02) 


1 -(S+O 9 ) a AV)2 v ( \ i 

-a 2± e V U / ^x 0 (-^ 02 , S ) 


(s € S 2i ) (52) 
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„ , v , , x c 2 a 2 ^lX^s) , . -2(s+a 2 )a/^ 


ik 1 


pi i 

\ E2±(^’ s )' 
J 0 


2 X 0 (-f,s) 

5 2 (s) a 2± e ( S+a 2) a A02 x 0 (-i/ 02 ,s) 


2 ( s+a 2) a /^ X Q (-p,s) (? 1 ; ,(m,s) 2 p dp 
S2 2 (M,s) £2 2 (p,s) c i cr l I/ 


02 


(O^si) (53) 


and 


1 , , . -(s+a-,)a/y ni 

|c 1 cr 1 ^ 01 J2 1 (i/ 01 ,s) a 1± e ^ lJ ' 01 


- J l±( y 0l) + 


k, 1 


2 X 0 (^ 01 ,s) 

« 2 (s) a 2± e ' (s+t ' 2)a/4 ' 02 x 0 (-^2,s)-^- 


f 1 e 2 ± (m,b) £ 2 ^ <72 ) a / ,1 x l: ,(-y,s) >t dM 

J 0 to£(n,s) fi 2^,s) (/* ■ ^Ol) 


’02 ~ y 01 


(s e Su) (54) 


The I m± and J m± terms in equations (51) to (54) which contain only integrations over 
the initial distribution are therefore known functions when f(x,p) is specified and are 
given by 

, v Ci a-. (s+a-iWt' . 

hd v ) = ^r F i±( -a ’^ s ) eV ; ^2^» s ) fi 2(^ s ) 

Li Li 


1 v^2My ( „ c \ 
2 k %Ks) X ° ’ 


Mil H' V 


§2(s) F 2± (a^ 02 ,s) e ( s+g 2) a Ao2 X Q (-i/ Q2 , S ) ^ ^ 


_ %(°°,s) 

+ <SRi) 


C 1 F 1± (-a,M,s) e( S+< ^ 2^ 
J 0 


02 

(p,s) 2p dp 


+ 6j(s) F 1± (-a,v 01 ,s) e ^ S+CT ^ a /^ 01 


X 0 (-p,s) c 2 a 2 i/ 


y 01 


x o<Moi> s ) y “ y 01 


(0^11) (55) 
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J 2±(^02) 


Il±(^) = + 
± 


+ 


and 


= ± 




al „ , , -(s+a^aAt , M d/x 

F 2± (a,M,») e <■ V / Xgi-w)—— 


| F 2±( a .‘’02’ s ) e"( S+a2 ) a / I '02 X 0 (-r 02 ,s) 


2 

%(°°,s) 


« 2 (°°,s) 


f 1 F u (-a,n,s) e ( S+<ri ) a / tl - JiSfe . 

Jf) X 0("^5 S ) (M " ^02) 


+ 6^(s) F 1± (-a,^ 0 i,s) e^ S+ai ^ a / I ' 0;L =r-r M \ 

1 1±v 01 ; x o(-^oi^ s ) (^01 “ ^02) 


(s e S 2i ) (56) 


-fs+cr-iWz' CoCTo , v -(s+a ? W 
F 1± (-a,^s) e ' l} ' - ^F 2± (a,^,s) e ^ 2 > / 


Si{(y,s) 


h 1 


£ F 2± (a,M,s) e ( s+a 2 ) a / M x q (-m,s) <p ly (M,s) 


_2 X 0 (-y,s) 

5 2 (s) F 2± (a^ 02 ,s) e - ( s+a 2) a / y 02 X 0 (-i/ 02 , s) 


02 


fi 2 (°°,s) 


f F 1± (-a,M,s) e ( s+tr l) a / , ‘ ■ 

J 0 1:t Xq(~h,s) (fi + 




5 l( s ) 


^01 

XO^OI’ 3 ) (*01 + *) 


(0 2 1 / < 1) (57) 


+ F l±(" a ^01. s ) l c l CT 1^01^l( y 0b s ) e 


(s+ai)a/^ 01 


1, 1 

-k: 


[2 X 0 (-^oi 5 s) 


_ / \ - ( s+cr 2) a /^ < ’ -u- / \ M du 

F 2 ±( a ,M,s) e ' 2; ' X 0 (-jz,s) -f- M 


M " ^01 


6 2 (s) F 2± (a^o 2 ,s) e ( s+g 2) a Ao2 X Q (-y 02 ,s) — - 


02 " "01 


%(°°,s) 


C 1 Fi (-a a s) e^ S+ai ) a / M - d/1 

J 0 X± ’ ’ Xo(-fi,s) (fi + y 01 ) 


« 2 (~ s ) 


(s e S n ) (58) 
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In these equations, the XQ m functions which were shown in reference 14 to be continuous 
across the curves C m in the s-plane (see appendix A and fig. 2) have been used. For 
two material media, one takes the ratio of these single -medium XQ m functions 


X 0 (z,s) 


Xq2(z,s) 

X 0 i(z,s) 


(59) 


where 


and 


x 0m( z > s ) 


{ v 0m - z ) X m( z > s ) 
(1 - z)X m (z,s) 


(s e S m jJ 
( s e s me) 


X m ( z j s ) 



(60) 


( 61 ) 


For Re(z) < 0, Xq(z,s) given by equation (59) is a nonvanishing analytic function of z 
and s provided s / ^-cr m , _cr m (l - c m )j, the branch cut of yo m (s), m = 1, 2. The 
quantity 


k = s^c^ - c 2 a 2 ) + ^ 03 (c x - c 2 ) (62) 

is related to the difference between medium 1 and medium 2 continuum solutions; several 
equivalent expressions for k are given in appendix C. 

Equations (50) to (58) determine the expansion coefficients A m± and a m± as 
follows. Recall that the inhomogeneous terms I m± and J m± are known functions for a 
given initial distribution f(x,/u). Equation (52) is used to eliminate a 2± from equa- 
tion (51) and upon using equation (50), one obtains an inhomogeneous Fredholm integral 
equation for the unknown coefficient A 2± . It is seen from equations (52), (53), and (54) 
that the remaining unknown coefficients are given in terms of A 2± and other known 
functions. However, one needs to know the analytic properties of the transformed solu- 
tion if/ ± in the s-plane in order to invert the Laplace transform. For part of this inves- 
tigation, another form of the solution is much more convenient. 

Complex Representation of i// ± (x,/i,s) 

In equations (50), the coefficient E m± (t',s) were introduced since they are the forms 
of the normal -mode expansion coefficients which are extendable to the complex plane. 


23 



(See ref. 21.) Thus, equations (51) to (58) can be written in a compact form valid for 
Re(s) > -a min = -min^rj,ff2). A brief outline of this extension to the complex plane is 
given in appendix E. 


E 2± (z,s) = I 2± (z,s) 


k ft 2 (°°,s) X Q (-z,s) 
c 2 a 2 n l( co J s ) 2ni 


l 


E 2± (z’,s) X 0 (-z',s) e 


-2(s+a 2 )a/z' 


^2( z 'j s ) (z' + z) 


dz' 


( 63 ) 


E 1± (z,s) =I 1± (z,s) ±^^E 2 ± (z,s) e 2 ( s+a 2) a / z 


c 2 CT 2 x o(“ z > s ) 2iri JC 


i Jr * 


E 2± (z’,s) X q (— z',s) e 2 ( s+p 2) a / z ^ 


« 2 (z',s) (z* - z) 


(64) 


! 2 ±(z,s) =^L 1± (-a,z,s) + 


k ^2 v / » 

2iri fij(<=o,s) 0 ’ 


r ^2± 

JpT CoC« 


L 2 ±( a > z ',s) x o(" z '’ s ) 


L Jc* c 2 a 2 fi 2 (z’,s) (z' + z) 


dz' 


%(°°,s) p L 1± (-a,z’,s) dz’ 


+ Sl 2 («>,s) J c . c 1 a 1 X 0 (-z , ,s) %(z’,s) (z’ - z) 


(65) 


and 


Il ± (z,s) = +L 1± (-a,z,s) e 2 ( s+a l) a / z ± £iJLL 2 ± (a,z,s) 


27riXQ(-z,s) 

n l(°°,s) 


L 2± (a,z’,s) X 0 (-z’,s) 
C’ c 2 ct 2 £2 2 (z’,s) ( z ' - z ) 

L 1± (-a,z',s) dz’ 


■S 

Jr t 


dz’ 


»l(°°,s) Li ± (-a,z’,s) dz’ 

+ ^2(°°> s ) Jc* c l CT l x o(" z, > s ) fi l( z ’» s ) (z* + z) 


The functions L m± (x,z,s) which appear in equations (65) and (66) are given by 


j m ±(x,z,s) = f 


i -(s+CT m )(x-X 0 )/z 


1 f ^ r / , du 

2 C m a mJ o f m ±( X 0’~^)]TT7 

2 C m a m | ^m±( x 0’^)^ _ z + z^m±( x 0> z ) ^m( z ’ s ) 


dx r 


(66) 


(67) 
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with 



and are analytic for Re(s) > -a m . These functions were introduced for z 
as 

, -(s+o rri ')x/i' 

L m±(x,^s) = F m± (x J y,s) O m (v,s) a m (v, s) e ' / 


( 68 ) 


v, 0^11, 
(69) 


in order to extend F m± to the complex plane. In equations (63) to (69), z does not lie 
outside the contour C' which encircles yQ m as shown in figure 3. The restriction 
Re(s) > -Ojrnn is discussed in the next section. 



That equations (63) to (66) reduce to equations (51) to (58) as all contours C 1 are 
collapsed onto the branch cut v e (0,1) due to f2 m (z,s) (refer to fig. 3) can be seen as 

follows. If s e S m ^, [p m (z)] ^ has a pole at z = i^Q m whose residue leads to a dis- 
crete term. When s e S me , £2 m (z) does not vanish. The continuum terms are simply 
those due to the integration around the branch cut. 
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The solutions )^ mc± (x,/j.,s) and \p m ^ ± (x,ii,s) can now be written simi- 
larly as 


<// 2c± (x,m,s) = 2^1 


IL 


E 2± (z',s) e 


' ( s+cr 2 )(a+x)/z' 


f 

JP! 


c ^2( z '»s) (z* - M) 

/ -(s+ao)(a-x)/z' 

E 2± (z’,s) e v 2J ' 


dz’ 


C* tt 2 (z',s) (z* + M) 


dz 


(l x | < a) (70) 


for Re(s) > ^min? 


^lc±( x ^,s) =/ 


/ . \ (s+CTi)(x+a)/z' 

_1_ C E 1± (z',s) e^ V d _ t 

27Ti Jqi fil(z',s) (z* + JU.) ~~ 
-(s+a 1 )(x-a)/z T 


±1 f E 1± (z',s) e 


n Jrr 


2tT1 Jqt ^(z t ,s) (z f - jLt) 


dz T 


IX < -a) 


(x > a) 


(71) 


for Re(s) > -o^: 


min> 


. , , 1 rp L 2± (x,z , ,s) r 

^2p±( x ’^’ s ) “ 27ri L n 2 ( 2 .,s) ( z « - ju) 62 ± J c , 


L 2 ±(-x,z',s) 

C* Q, 2 (z’,s) (z' + ju) 


dz 1 


(lx | < a) (72) 


for Re(s) > -a 2 , and 


27Ti 


Ll ± (x,z',s) 
c , %(z',s) (z* - jut) 


f 

Jr* 


dz' 


L 


M ± (x,z’,s) ± Li ± (-a,z',s) e 


-(s+o-i)(a- x )/z' 


%(z',s) (z* + ju) 


dz 


(x < -a) 


i//l p± (x,/r,s) =/ 


(73) 


2vi 


-(s+a,')(a+x)/z' 

C L l±( _a ’Z’ > s) e v ' t M.(x,z',s) 

\ — ^ _ , . > , . 5 dz’ 

Jr* 


L 


C* «l(z',s) (z* - ju) 

L 1± (-x,z',s) 


C* %(z',s) (z’ + JLt) 


dz 1 


(x > a) 


for Re(s) > -cr-^. 
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The functions M ± (x,z,s) are also integrations over the initial distribution fj ± (x,ju.) 
and are given by 


M 


p~cl 

.(x,z,s) = - l 


(s+a 1 )(x+x o yz 


f c i a i C f i 


and 


M. 


‘i Ciai I 0 f l±( X 0’" M )/^ + z f l±( X 0’" Z ) fil(z ’ s) 

;(XA s) . j; £ , uM ^ 

"| c l a l J f l±( x 0’-^)/^ + i f l±( x 0^ z ) %( z > s ) 


dxQ (x > a) (74) 


dxg (x < -a) (75) 


for Re(s) > - 0 } and z not outside C'. Again, the discrete and continuum terms which 
appear in equations (32), (39), (44), and (45) are due to the zeros and branch cuts of 
J2 m (z,s) which appear in the integrands of equations (70) to (73). 


PROPERTIES OF TRANSFORMED SOLUTION 


General Properties in s -Plane 

Analytic properties of ^(x,^, s) as a function of s must be investigated before 
the time -dependent solution ^(x^t) can be recovered according to the inverse Laplace 
transformation given by equation (7). To do this, the behavior of \j/ ± in some right-half 
s -plane is required. Before looking at the details, a review of some results of earlier 
cited work in which Case's method was used is in order. In these works, the analytic 
properties of the functions of s such as z/Q m , anc * the various X-functions are 

given. 

In the semi-infinite medium problems considered in references 14 to 16, expan- 
sion coefficients could be found explicitly and this fact aided in the extraction of the 
s -dependence of the transformed solutions. These solutions were found to contain the 
branch cuts of !^g m (s) so that the integration contour of the inverse Laplace transfor- 
mation was deformed around these branch cuts. For the slab problem solved in refer- 
ences 12 and 13, expansion coefficients could not be found explicitly but the theorems of 
references 6 and 7 gave the analytic properties of the transformed solution in the s-plane. 
In that problem, the solution does not contain the branch cut of i'q(s) even though vq(s) 
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appears explicitly in it. Instead, there are a finite number of poles at values of s, for 
example, Sq, . . ., s N , which lie on the branch cut of i'g(s), that is, on the real s-axis. 
These poles contribute a sum of residues as the integration contour is moved to the left 
of them in the s -plane. Furthermore, in these previously solved time -dependent prob- 
lems, there is a real number, for example y^, such that the integration contour cannot 
be deformed into the region Re(s) < y^ for arbitrary values of x. The present trans- 
formed solution should exhibit similar properties; that is, xj/ ± may not be analytic for 
Re(s) less than some number y\ when x is arbitrary whereas for Re(s) greater 
than yj, it should be analytic except for poles and/or branch cuts. Such singularities 
probably occur where ^Q m (s) has its branch cut. 

First note that for arbitrary initial distributions f(x,p), i// ± (x,p,s) is not analytic 
for Re(s) < -cr m j n . This statement is true since each of the inhomogeneous terms I m± 
of equations (63) and (64) contains both L 1± and L 2± as can be seen from equa- 
tions (65) and (66) and one of the two will not be analytic for Re(s) < -a min = -min 
In particular, note that for |x| > a, \pi ± (x,n, s) never appears to be analytic for 
Re(s) < -a m j n . However, for special cases of material properties and initial distribu- 
tions, ^ 2 ±( X )M, S ) can be shown to be analytic for -02 < Re(s) < -oq except perhaps for 
poles. 

Consider now the behavior of for Re(s) > -Ojnin* Note that the transform 
plane for the present problem must be taken as a superposition of two "single-medium" 
planes, that is, one for each material medium in the problem. The expressions (32), (39), 
(44), and (45) for the transformed solution were not defined for s e C m and outwardly 
appear to be discontinuous at s e C m . However, this is not the case. The complex 
representation of E m± given by equations (63) and (64) shows that such coefficients are 
continuous across the curves C m . Thus, it is seen from the representation of xf/ ± 
given in equations (70) to (73) that i// ± is indeed continuous across the curves C m . 

The Associated Eigenvalue Problem 

It is convenient to introduce at this time the solution of the associated eigenvalue 
problem, that is, the solution of equation (18) subject to the boundary conditions (19) and 
(20) with f m± (x,/i) s 0. Such solutions, denoted with a bar, have the form 
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\p ± (x,i±,s) =< 


b l±^_j ; 01 (x,M,s) fii(s) + B 1± (-y) ^l(_„)(x,JLl,s) dy (x < -a) 

[% 2 (X ’ M,S) ± ^-^ 02 ^ x,m ’ s) 5 2 ( s ) 

+ J o B 2± (y)[^2y(x,M,s) ± ^ 2 (_ y )(x,/i,s)] dn (Ixl < a) 

_ pi 

±b l±>/V 01 (x,M,s) 6 l( s ) * J B l±(-v) dv (x > a) 


(76) 


where obviously B m± and bi ± can be obtained from the E m± given by equations (63) 
and (64) in the case f m± (x,/i.) = 0. It will be seen later that the solution x{/ ± has poles 
at those values of s for which the associated eigenvalue problem has nontrivial solu- 
tions. In appendix F, it is shown that as the slab thickness becomes very large, this 
eigenvalue problem has only trivial solutions for Re(s) > -cr 2 except perhaps on the 
branch cuts of ^Om^' When the slab thickness is not large, one still expects that if the 
eigenvalue problem has nontrivial solutions for Re(s) > -ct 2 , they occur only when s is 
real. This statement has been proved rigorously by using the approach of Lehner and 
Wing (refs. 6 and 7) for several problems which are special cases of the present problem: 
the bare slab considered in references 6 and 7 and the slab surrounded by pure absorbers 
considered in references 8 and 9. In all these problems, there is no scattering in the 
reflector and, therefore, no branch cut of ^(s). As already indicated, the Xq(z,s) 
function contains the branch cuts due to both z-'oi(s) and fQ 2 (s) and these branch cuts 
lie on the real s-axis from -a m to -cr m (l - c m ) and may or may not overlap depending 
on the values of material properties. Note that has been taken less than unity and 
this insures that the branch cut of vqi lies entirely to the left of s = 0. In previously 
solved time-dependent problems, singularities of the transformed solution always occur 
where the ^m( s ) has branch cuts. Since the analysis of appendix F indicates that for 
large values of the slab half -thickness a, the singularities of \ p ± for Re(s) > -o^n 
also occur where the ^ 0 m( s ) have branch cuts, it will be assumed for all values of a 
that the singularities of i{/ ± occur on the branch cuts of t'o m (s). In any case, it is 
shown that the only other singularities of i// ± , Re(s) > -a m ^ n which could occur off the 
branch cuts of i^o m ( s ) are poles, whose residue could readily be added to the time- 
dependent solution. 
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In order to determine the behavior of xf/ ± on the branch cuts of ^o m ( s ) one f irst 
considers \j/ ± in the region s e S-q D S 2 p For this region, the expansion coefficients 
are given by the following equations. (See appendix G for a brief description of the man- 
ner in which these equations were obtained.) 


b 2 ± (m) = 


k ^2K s ) 1 

( y 02 " l±2 ] 

1 h 2 (jLt) 

h 2(^02) x h 2(“^02) 

± 2 %(°°,s) | 

W - » 2 ) 

1 n 2 (m) 

_M + ^02 M “ ^02 


+ y o 1 B 2±M h 2M? fL. 


(0 ^ jLl 2 1) (77) 


B l±(“^) +^ B 2±(M) e^ 1 


k s ) h l(^) 


•2 S ) N iC/x) 


± 52(^02} + f 1 Bn Av) W,) d- 

V--V 02 li + ^02 Jq 2± 2 c l a l^ 


(78) 


and 


+ h l(-^0l) b l± = h 2( I/ 02) ± h 2("^02) + ( u 02 ~ v 0l)§ 0 B 2±M h 2^ ~ ^^ ~2 (79) 


where 


h 2 (w) = 


X 2 (-oj,s) -(s+a 2 Wu> 

\ e 11 

Xj_(-w,s) 




hi (a>) = 


f2l(°°,s) Xi(-w,s) (s+cr^a/w 
« 2 (°°, s) a, X 2 (-w,s) 6 




N m (M) = s ) fi m(^ s ) 

In addition, the eigenvalue condition 


0 _ h 2 (^02) 

^01 + ^02 


h 2(- y 02) 

^01 " v 02 


lo B a± w ■’aMTTToI 


(80) 


(81) 


must be satisfied. As noted in appendix G, this equation is an additional constraint on the 
solutions of equation (77). Since the eigenvalue condition (81) has different limiting 
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values as s approaches the branch cut of ^(s), it is concluded that there are only 
trivial solutions of the associated eigenvalue problem on the i'gi(s) cut. When s 
belongs to the branch cut of ^(s) which is not also part of the i'qi(s) cut, that is, 
when He(^o2) = Im (^oi) = it appears that nontrivial solutions of the associated eigen- 
value problem may exist. From Bowden's results for the bare slab (refs. 12 and 13), it 

is expected that equations (77) and (81) are satisfied only at isolated points ^s n ^ . In the 

limit C202a — these points lie on the branch cut of ^(s), that is, the s n are real. 
The "thick-slab" eigenvalue condition is seen from equations (77) and (81) to be equa- 
tion (81) with B 2 ± (m) = 0. 

If material properties are such that -c^ < -0|, then part of the branch cut of 
fQ 2 (s) lies in s e S 2 i fl S^ e . In this region however, s < -cr m i n = -o^ and for such 

values, the solution i// ± (x,/i.,s), |x| > a, that is i^ ± , is not bounded as |x| — °°. How- 
ever, i// 2 ± may have nontrivial solutions on such a part of the branch cut of 1^2 (s). 

The equation for B 2 ± and the additional constraint for this region are (again, see 
appendix G for some discussion) 

r („\ - . k ^2(°°> s ) ( 2 ,.2\ h 2^ X l(-M,s) | h 2 (^Q 2 ) x l(~^Q2> s ) 

2± 2 %(°°,s) \°2 /N 2 (m) X 01 (-/i,s) M + ^02 x 0 l( _i; 025 s ) 


h 2 (~ ^02) X 1 (^02 j s ) 

H - v 02 x 01 ( I, 02,s) 



x l(~^ s ) dv 
X 0 l(-^,s)^+M 


(0 = jla = 1) (82) 


and 


0 = h 2(1,02 * h2( "' 02) ^fe| + £ B2±M h2 °' ) 


Xi(-^,s) 

Xqi(-^s) 


dv 


(83) 


It will be seen later that the zeros of equation (83) can, under some conditions, be poles 
of ^2± and therefore may contribute discrete modes in ^(x,/Li,t), |x| < a. For this 
reason one is interested in where these zeros lie. They will be referred to as 
pseudo-eigenvalues. 


Relationship Between i// ± and if/ ± 

It is now shown how the solution of the associated eigenvalue problem i// ± is con- 
tained in the inhomogeneous solution if/, by following a procedure similar to that of 
Bowden and Williams (ref. 13). In appendix H, it is shown that the original expansion 
coefficients of equations (32) and (39) can be written as 



and 


a l± = 


a 2± + | F 2±( a ^02, s ) 


b l± + b l± 


( s e S li n S 2i) 


where B m± and b^ ± are given by equations (77) to (79). The coefficients B m± 
b^ ± are given by 




k n 2 («,s) (» M » - r 2 ) h 2 (K) \A 

2 OjK s) ^2 . Jj N 2 (k)SJ 0 2 = 0) 2<m) U + 1' 


h 2 (^02) h 2 (~^02) 


+ 2 F 2±( a ’ y 02’ s )|y + uq 2 n v ~ ^02 
+ f F 1± (-a,jLi,s) h 1 (/x) 


+ J F 2 ± (a,M,s) 


n 

"01 - V - 2 2 ?> 2 j ,((i,s) 


d /lx 


k ^2(°°>s) h l( y ) 


" ± 2 «i(“,s) Ni(v) \2 F 2±( a ’^02) s ) 


h 2 (^02) _ h 2 (~ ^02) 


v ~ ^02 * + v 02 


pi r -1 2( Pi 7 ,(MjS) 

J n L B 2±^) + F 2± (a,^,s)j h 2 (m) ■■ e ' glt , d/x 


P 1 Kf 

T \ F 1 ± (-a,M,s) hj(M) 7 — 5 - 

J ° ("02 


("oi 2 -^ 2 ) dn 


M 


2 _ „2^M + * 


F 1 ± (-a^,s) -^^F 2± (a.,v,s) e ^ 1 °^ V 


and 


+ h l(- I, 0l) b l± - F ±(- a ^0b s ) 




(84b) 

and 


(85) 


(86) 


(87) 
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The coefficient a 2 ± + ^ ^ r 2±( a ’ t; 02) s ) i s given by 
„ , Ip _ ~ v 0 lPl ± + ^ 2 ± 

L 2± 2 2 « ’ 02 ’ il ' “2±) 

In these equations, a! 1± , a! 2± , Pi ± , and /S 2± are 

a l ± = h 2 (^ 02 ) * h 2 ("^ 02 ) + (^02 " ^ 0 l ) b 2 ±^) h 2 ^) ~ 9 ^ 

^0 n* - v. 


* ~ V 01 


a 2± = y 02 h 2 (^ 02 ) + y 02 h 2(-^02) + ( y 0 2 “ V 0l) ^ B 2± -M h 2 (n) 

01± = f F 2±( a ^02’ s )[^2( z/ 02) + h 2( -I/ 02)J 

+ (^02 - ‘'of) J 0 [ B 2±(^ + F 2 ±(a,M,s)] h 2 ( M ) 2 

M ~ ^oi 

* I 1 " 01 “ V oi) C, F l±(- a .*M^ s ) h l(M) 2 2 


li dju. 

71 T~2 

M ^oi 


^ 


± |Fi ± (-a^ 0 i,s) hi^Qi) + F 1± (-a,-^ 01 ,s) hi(-^ 01 2 
^2± = | F 2±( a > i/ 02’ s ) ^02 {^2 (^ 02 ) ± h 2( _I/ 02)] 


+ (^02 2 - ^OlK ^2*^ + F 2±< a ’^ s >l h 2^) 2 

u ^ - ^01 

*(*01 - v 02 2 ) t F l±(- a 5 M,s) hi(M) -~^~2 


li 2 - v Q l 


’01 F l±(“ a ^01’ s ) h l( v 0l) ~ ^01 F l±( _a »" y 01» s ) h l(“^0l)] 


33 



In terms of these quantities, the solutions i// m± can be written as 


i// 2± (x,ja,s) =ja 2 ± + |F 2 ± (a,v 02> s) ^ 2± (x,ii,s) + j* B 2± (u) fy 2l> (x,n,s) ± xp 2 ^(x,ii, sj] d? 

+ £ l t / 2v( x ’L L ’ s ') ± f 2±(- x >^ s ) ^2(-i^)( x ’^> s il dz/ 

+ ||f2±( x ^ 02’ s ) ± F 2±(“ x > -I/ 02» s )] % 2 ( x »M,s) 


+ =•= F 2±( _x > ^02 j s )J *-v q2 ^,s) 


( I x I < a) (91) 


and 


xly 1± (x,ii,s) = ja 2± + |- F 2±( a > I/ 02j s ) <^ 1± (x,M,s) 

± [t>i ± - F ± (-a,i/ 0 i,s) + F l±(- x »-^0b s )] ^i^ 01 ( x ^,s) 

± Fi ± (- x ^ 0 l,s) ris_ VQ1 (x,iJ.,s) 

± [B 1± (-v) - F ± (-a,^,s) + F 1± (-x,-^,s)| i // lt ,(x,M,s) 

±j* F i±(- X ,^ s ) ^l(-z,)( x >M> s ) dv 


(x > a) (92) 


The solution i//j ± (x,|Li.,s) for x < -a has a similar form. In these equations, 
^ m ±( x )M>s) are the parts of i// ± (x,/x,s) which are given by equation (76). Equa- 
tion (79) is written in terms of a-^ ± as 

+ h l(-^0l) E l± = a l± 


(93) 
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Properties of i p ± on the Branch Cuts of ^o m (s) 

Consider now what happens on the branch cut of i'qj(s) where vq± = i | l| ^ or 
Im(s) = 0” and vq± = -i^oil f° r Im(s) = 0 + . From these equations, it can be shown 
that the quantities B 2 ± , B^ ± , B 2 ± , Bi ± , ai ± , ot 2 ± , Pi ± , and @ 2 ± do not contain 

the branch cut of rgi(s). Equations (93) and (87) show that bi ± and bi ± have branch 

cuts due to that of t'gi(s). Equation (88) indicates that a 2 ± + ^ F 2±( a,t, 02> s ) has the 

branch cut due to i'qj.(s) unless &\ ± /<X2± is e( l ua l to $l±/fi2±' 111 general, this state- 
ment will not be true since Pi±/@2± depends on the arbitrary initial distribution f ± ( x , m ) 
whereas &\±/ a 2±. does not - Therefore, it is concluded that both ip^ ± and i// 2± con- 
tain the branch cut of z'oi(s). 

On the branch cut of ^(s), the <l ua ntities B 2± , Bj ± , bj ± , /3j ± , and 0 2± are 
single valued. Since the quantities a;j ± and Q! 2± of equation (89) are related above 
and below the branch cut of ^Q2( s ) by 



it follows then from equation (88) that on that part of the branch cut of ^q 2 (s) which is 
not also part of the rg^(s) cu ^ that * s > tor ■^ e ( l/ 02) = Im^oi) = one has 


a 


a 2± + l F 2±( a -’ v 02> 



(95) 


if the denominator on the right-hand side of equation (88) does not vanish. It is seen from 
equations (76) to (79) that for this same region, 



Hence, the product 


(96) 


a 2± + | F 2±( a ^02 J s ) 


4/ {x, ii, s) 


(97) 


which appears in if/ does not contain the branch cut of ^q 2 (s). However, the denomi- 
nator of a 2 ± + ^ F 2 ±( a , ^02?®), namely " a 2±> equivalent to the eigenvalue con- 

dition (eq. (81)). Thus, if the associated eigenvalue problem has a nontrivial solution at 
s = s n , Re(s) > -CT min then, if/ has a pole there. 
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The analytic properties of the transformed solution ij/ ± (x,iA,s) may be summarized 
as follows. For arbitrary initial distributions f ± (x,p), i f/ ± is not analytic to the left of 
Re(s) = ■o'min in the s-plane, whereas to the right of Re(s) = -cr min it is analytic except 

for the branch cut along -0i(l - (due to the branch cut of Z'oi(s)) it 

CT m in > crj(l - cj) and poles at the values of s at which the associated eigenvalue prob- 
lem has nontrivial, solutions ij/ ± . It has been assumed that for arbitrary slab thick- 
nesses a, these poles, if they exist, lie on the branch cut of vq2 ( s ) since this result is 
the rigorous one obtained by others for several special cases of the present problem and 
obtained herein for the case when C2<J%a. is large. For special values of material prop- 
erties and initial data, i// ± (x,/a,s) for |x| < a (that is, j/^ 2±) ma y he analytic in the 
region -02 < Re(s) < -Oj except perhaps for poles. 

RECOVERY OF TIME -DEPENDENT SOLUTION 

The time -dependent solution 'Zr(x,ii,t) is obtained from the inverse Laplace 
transformation (eq. (7)) where y is to the right of all singularities of i//(x,p, s) 
in the s-plane. From the preceding analysis, it is expected that one can choose any 

y > max^-aj(l - cj), -03(1 - C£)^. In order to show the time dependence of the solution 
4>(x,p,t) more explicitly, the inversion contour should be deformed as far as possible to 
the left in the s-plane by making use of the analytic properties of i//(x,p,s) obtained in 
the previous section. 


Behavior of 1 f/ ± on the Contour Re(s) = y 

The behavior of rp ± on the contour Re(s) = y must be examined as |s| — 
(See fig. 4.) This contour crosses both of the curves C-^ and C2 and it has been 
shown that i p ± is continuous across these curves. As |s| — °° on such a contour, 
s e S le n S2 e and it is shown in appendix I that i// ± behaves as follows: 




f 2±( x 0^) + °(i 


dxr 


+ 





( 98 ) 


for |x| < a and p > 0, 
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^2: 


i<x,m ’ s) - r r e 


-(s+a 2 ) 


x 0“ x 

iMl 


f 2±( x 0^“ l^-l) + °(|) 


dxQ 


a-x 


. ( Cr 2 _Cr l J /\|/i.| I - 


Im| 


I 



0" X U 


m 


fl ± (x 0 ,-|Ml) + °(|) 


dx r 


for x| < a and pi c 0, 


( a 2" a l)( ix ) pa -(s+02)(— f \ /i\ 

^l ± (x,M,s) £ e V d ) |f 2± (x 0 ,|i) + °( 5 ) 


Mr e -N^ 

^ J - 00 


[ 1± 


(x 0 ,/x) + o|j 


dx r 


dx r 


a — 

for x > a and pi > 0, and 

>.M - ]M” e "( S+<7 l)(w) t 1± (x 0 , -|fi|) + 00) 


x-: 


x 0' 


'l± v 


dxr 


(99) 


(100) 


(101) 


for x > a and pi < 0. Expressions similar to equations (100) and (101) are obtained 
for x < -a. It is seen that i/ ± is not necessarily 0^. However, the parts which are 
not can be easily inverted as follows. Define for all s the function i// u± (x,pi,s) as that 
part of each of equations (98) to (100) which is not o(^. It is shown in appendix I that 
upon making the substitution 


X - Xq = pit 

xq - x = |p|t 

i// u± (x,pi,s) can be written as 


(pi > 0) 1 

(pi < 0) f 


(102) 


^ u± (x,M, 



e -st 



(103) 
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that is, the parts of \p ± which do not behave as 0[±J as |s| — Re(s) = y can be 
inverted by inspection. The solution ^ u± (x,iU,,t) is given by 


f 2± (x-fit,jLi) 


* u± (x,/i,t) = <j 


- CT lt -( a 2- a l)(^) 


e e 


fl ± (x-/it,/i) 


(‘ <a r 


t > 




(104) 


for |x| < a and /i > 0; 


T u± (x,/i,t) =/ 


for |x| < a and ju. < 0; 


e ° 2 * 


f 1± (x-/xt,/x) 


t < 


t > 


a - x 

H 


a - x\ 

W / 


(105) 


^ u± (x,M,t) =< 


r -a-ft 

e 1 f 1± (x-Mt,M) 

J.-v 


f 2± (x-/it,M) 


-^t -( a 2- a l)^f). , . , 

e 1 e f l±( x -Mt,|U.) 


(<^) 


x-a < j ;< x + a 






-) (106) 


t > 


x + a\ 

M J 


for x > a and p. > 0 and 


-CTlt 

* u± (x,M,t) = e f 1± (x-Mt,/i) 


(107) 


for x > a and p. < 0. That \& u± describes the motion of uncollided neutrons from the 
initial distribution can be seen by direct substitution; that is, ^ u± satisfies the equation 


8* 8* 

_p +M _H± + c,( x )* u± = 0 


(108) 


In the limit t — 0, note that 


* u± (x,/li, 0) = f ± (x,p) 


(109) 
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For arbitrary f(x,/i) which vanishes as |x| — °°, !// u± (x,jn,s) given by equa- 
tions (103) and (104) to (107) is an analytic function of s for Re(s) > -cr min for almost 
all x and i±. If fi ± s 0 (f 2 ± = 0), then i f/ u± is an analytic function of s for 
Re(s) > -o^ (Re(s) > -cr-jJ. Therefore, the function 4> ± (x,/i.,s) defined as 


4> ± (x,m,s) h i// + (x,/i,s) - i// u± (x,/i, s) 


(Re(s) > (110) 


has the same analytic properties as xl/ ± in the right-half plane Re(s) > -cr m i n except 
that it is og) as |s| - If i// ± has a branch cut along ( _cr m inj _a l(l “ c l))j that is, 
if a min > ai(l - c x ), then 


on the branch cut. Similarly, if i f/ has a pole at s = s n , then 


( 111 ) 


Residue 


(<f> ) = Residuef 

\ ±/s n 


( 112 ) 


Explicit Form of ^ ± (x,/i,t) 

The definite parity parts of the time-dependent solution therefore can be written 
from equation (7) as 


i py+i 

tf ± (x,#i,t) = * u± (x,M,t) + \ 


2S) r -l„ eS ‘ dS 


(113) 


By using the analytic properties, one can deform the contour to the left and obtain, in 
general, 


U±(x,M, t) + 

I 



3-S n 

1 

f^iC 1- 

■ C l)j 

27ri 

J-a . 

i 


min 


1 

f“ a min 

+ ioo 

27ri 

J_ a . . 

u mm 

-ico 

1 

lim f 

lb 

27ri 

P -0 J c 

, Y ± 
’ r\ 


e s ^ ds 


+ 27 tI 0 ™O Ip ^± (x ’ m ’ s) eSt ds ( _cr min < -CT l( 4 “ c l) < s n) ( 114 ) 
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where Cp is a small circular contour of radius p with center at s = -o^(l - c^). All 
the contours are indicated in figure 4. Generally, the point s = -o^(l - c^) will not 
satisfy the eigenvalue conditions (eq. (81)) and the contribution from the contour Cp 
vanishes as p — 0. If, however, s = -crj(l - c-jj happens to satisfy equation (81), the 
contribution from the contour Cp has the form of a discrete residue term. Details 
concerning this point are discussed in appendixes I and K. 


Deformed Original 



Figure 4.- Integration contours for the inverse Laplace transformation of 
y(x,p,s) when vq^(s) branch cut is embedded in the vq 2 (s) branch 

cut, C 2 > 1 , and the initial distribution is arbitrary. 

Equation (114) is the solution of the time -dependent problem written in a form in 
which the uncollided portion of the initial distribution f(x,p) has been separated. For 
arbitrary f(x,p) the contour cannot be deformed further to the left. It is indicated 
in the final section that this solution reduces to those obtained previously by others for 
special cases of the present problem. 

This section is concluded by indicating the form of some parts of equation (114). 
The uncollided term 4i r u± (x,fi,t) is given explicitly by equations (104) to (107). The 
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form of !// ± (x,/n,s) on the branch cut -crj^l - c^)j was given in the previous 

section. From those results, it follows that on this branch cut, |V ± (x,ju.,s)] - (j^ ± (x,/i,s)] 
can be written from equations (91) and (92) as 

[>M x ’ M ’ s) ] - [^ ± (x,/i,s)| = jja 2± + |F 2± (a,^ 02 ,s)j - |a 2± + | F 2± ( a ,^02’ s )] J^ 2 ± (x ^> s) 

(115) 

for |xl < a and 


[^±(x,M,s)] - [^ ± (x,m,s)J = a 2± + |F 2± (a,^ 02 ,s)] ^ 1 ± (x,m,s)] 

- [ a 2± + | F 2±( a > ^02 ,s )] [^i ± (x,/i,s)| 

± f[ b l± " F ±(- a ^01’ s 3 %/ x ^’ s ) 

- b i± - F ±(- a ^oi» s j] ^ 01 i^ 


(116) 


for x > a, where ja 2± + | F 2±( a > l '02> s )| is g iven b y equation (88), ^ m ±(x,M,s) by 

equation (76) and jbj ± - F ± (-a,i'Qi,s3 by equation (87). The solution i// ± (x,/i,s) has 

poles at s = Sq, . . ., s N because of the poles of ja 2± + F 2±( a > v 02> s ^*'o2 + ^ / ^' A §ain, 

from the results given in the previous section, it follows that 


Residue |^/ ± (x,/i.,s) e s ^ 


Sn 


= e Snt L ± (x,ii,s) 4 1+1)/2 _ ) Residue<j4 +1)/2 [a2 ± + f F 2±( a ^02> s ) 




% 




(117) 


Js n 


Note that the factor u q 2 + ^^ is introduced so that ^± V Q2^ + ^ ^ an< ^ 

|*2± + | F 2±( a ^02 4 4 +1 )/ 2 are single valued on the branch cut of ^q 2 (see eqs. (95) 

Snt 

and (96)). These terms have an exponential time dependence e and the implicit 
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equations, namely, equations (77) and (81), from which the eigenvalues <(s n } can be 
computed have been obtained. Since information concerning the behavior of eigenvalues 
(that is, number, location, etc.) as a function of material properties is not readily 
obtained analytically from such expressions, a numerical study of real time eigenvalues 
has been made and the results are discussed in the next section. 

CALCULATION OF TIME EIGENVALUES 


First note that the eigenvalues and pseudo -eigenvalues depend on five parameters 
( c l> a l> c 2’ a 2> anc * a ) and therefore many numerical computations would be required 
in order to determine the specific dependence on each parameter. It will be seen that the 
bare-slab results of references 12 and 13, the theorems of reference 9 for slabs sur- 
rounded by purely absorbing media, and some observations of the present numerical 
results for a few reflected slab cases allow some conclusions about the behavior of 
eigenvalues for reflected slabs as a function of the slab half-thickness a to be drawn. 
However, rather than compute eigenvalues { s n} in terms of cj, oq, c 2j <? 2 , a nd a > 
one defines a nondimensional variable £ and nondimensional parameters ctr, cfq, and 
A as 


C = 


S + cr 2 
c 2°2 


c l g l 
R c 2 a 2 


ctt, = 


°D = 


g l " g 2 
c 2 g 2 


A = c 2 a 2 a 


j 


(118) 


In terms of these quantities, the branch cut of t>Q 2 becomes the real interval (0,1) and 
the branch cut of uqi becomes the real interval (-ctd, - g D + g R)’ Since a m and c m 
are nonnegative, it follows that 

-»d s 4 (U9) 

where the equality holds only if = 0. Also c^ has been restricted to less than unity 
so that -op) + S 1 implies that c 2 < 1. Obviously, = 0 when the reflector is a 
purely absorbing medium or a vacuum and = 0 when the total macroscopic cross 
sections of the two media are the same. It has been shown in the previous section that, 
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in general, the inversion contour can be deformed to the left only as far as Re(s) = -cr mi 
which corresponds to Re(£) = max^-aj),0). However, there are no eigenvalues on the 
branch cut of z'qi so the region of the real £-axis where the eigenvalues <^ n J> should 
appear is 


max (-or D + a R ,o) < ? n < 1 


( 120 ) 


This interval corresponds to s e n S 2 i and equations (77) and (81), written in terms 
of the quantities of equations (118), are solved numerically to obtain the real eigen- 
values for specified a R , ct-q, and A. The pseudo-eigenvalues are obtained 

numerically by solving equations (82) and (83) also written in terms of the quantities of 
equations (118). In addition, numerical results are obtained in the thick-slab approxima- 
tion, that is, equation (81) with B 2 ± (iu) = 0. Details concerning numerical procedures 
and computational equations are given in appendix J. 

The time dependence of discrete modes is seen from equations (114) and (117) to be 


s_t ( c 2^n -1 ) a 2 t 
e n = e 


( 121 ) 


Now £ n = -a D + a R implies that s n = -cr^l - CjJ ^ 0 since c^ < 1 and the equality 
holds only if Oj = 0. Therefore such £ n values correspond to time -decaying modes 
regardless of the value of C 2 - For values of £ n within the interval (120), the time 
decay or growth depends on whether C2? n is ^ ess than or greater than unity as can be 
seen from equation (121). A discrete mode represents a critical system if C 2 ? n = 1. 

The largest eigenvalue £q with an even parity eigenfunction corresponds to a critical 
slab problem with parameters 

c slab = ^ 

°R 

c reflector ~ C 0 + a D / (122) 


°slab a critical “ 

) 

where a- cri ^ ca j is the critical slab half -thickness. It is generally known (see, for 
example, ref. 2) that the critical radius for a bare sphere = 0) can be obtained from 
the largest slab eigenvalue with an odd parity eigenfunction. That is, when is 
used in equations (122) in place of £ 0 , the a cr ^i ca l is the bare-sphere critical radius. 
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Comparison With Published Numerical Results 

Many different combinations of material parameters could be considered, but here 
the study of the eigenvalue behavior is restricted to the case of overlapping branch cuts. 
As cr^ departs from zero, one would like to see how the eigenvalues depart from those 
previously reported (refs. 12 and 13) for a bare slab. A comparison of the present 
eigenvalues <^ n j> for vacuum reflectors, that is, cr R = 0, with those of reference 12 is 
given in tables I and II. Results generally agree to three figures for slab half- 
thicknesses A from 0.4 to 20. In table II, eigenvalues calculated in the thick-slab 
approximation are also shown for bare slabs. For slabs with half -thicknesses A > 1, 
the thick-slab approximation generally agrees with the numerical solution of the exact 
eigenvalue condition to three figures. This agreement can be seen from table III where 
such results are compared as departs from zero with Op = 0. From the bare- 
slab results = o) of tables I, II, and IE, critical slab half-thicknesses are obtained 
from £q by using equations (122). These values are compared with the critical slab 
half-thickness results of Mitsis (taken from ref. 2) in figure 5 (open symbols). Closed 
symbols give critical sphere quarter -diameters obtained from equations (122) and ^ 
whereas Mitsis' critical sphere results are taken from reference 22. The agreement is 
good to the scale of the figure. For or = 0, the eigenvalues £q and have also 
been compared directly with numerical bounds computed by Mullikin (ref. 23) for bare 
slabs and spheres and again the agreement is good. Critical half-thicknesses for slabs 


TABLE I.- EIGENVALUES (? n ) FOR BARE SLABS 


n 

Eigenvalues for - 

A = 1 

A = 5 

> 

11 

o 

A = 15 

A = 20 

Present 

Reference 12 

Present 

Reference 12 

Present 

Reference 12 

Present 

Reference 12 

Present 

Reference 12 

0 

0.703 

0.705 

0.975 

0.975 

0.993 

0.993 

0.997 

0.997 

0.998 

0.998 

1 



.897 

.897 

.971 

.971 

.987 

.987 

.992 

.992 

2 



.762 

.762 

.935 

.935 

.970 

.970 

.983 

.983 

3 



.560 

.560 

.883 

.883 

.946 

.946 

.969 

.969 

4 



.276 

.276 

.814 

.816 

.915 

.915 

.952 

.951 

5 





.728 

.727 

.877 

.877 

.930 

.930 

6 





.621 

.621 

.831 

.831 

.904 

.905 

7 





.493 

.493 

.777 

.777 

.874 

.874 

8 





.340 

.340 

.714 

.714 

.840 

.840 

9 





.157 

.157 

.642 

.644 

.800 

.800 

10 







.560 

.560 

.756 

.756 

11 







.467 

.467 

.707 

.707 

12 







.362 

.361 

.653 ! 

.653 

13 







.243 

.243 

.593 , 

.593 

14 







.110 

.110 

.526 ! 

.526 

15 









.453 

.453 

16 









.373 

.373 

17 









.286 

.286 

18 









.190 

.190 

19 




L . 





.084 

.084 
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TABLE II.- EIGENVALUE FOR THIN BARE SLABS 


Slab 

thickness, 

A 

Eigenvalues determined by - 

Thick-slab 
approximation 
(eq. (J19)) 

Present 
(eq. (J6)) 

Bowden 
(ref. 12) 

1.0 

0.702 

0.703 

0.705 

.8 

.612 

.615 

.615 

.6 

.473 

.483 

.483 

.4 

.244 

.282 

.282 

.2 

(*) 

.043 

.048 


*No solution found for £ > 0.001. 


TABLE HI.- EIGENVALUES C 0 AND <1 

[°D = °] 


FOR THIN REFLECTED SLABS 



Eigenvalues for - 


A = 0.4 

A = 0.7 

A = 1.0 

a R 

Co 

Co 

Co 


Thick-slab 

Present 

Thick-slab 

Present 

Thick-slab 

Present 


approximation 

solution 

approximation 

solution 

approximation 

solution 


(eq. (J19)) 

(eq. (J6)) 

(eq. (J19)) 

(eq. (J6)) 

(eq. (J19)) 

(eq. (J6)) 

0 

0.244 

0.282 

0.550 

0.556 

0.702 

0.703 

.2 

.376 

.398 

.600 

.604 

.727 

.728 

.4 

.512 

.522 

.661 

.663 

.759 

.759 

.6 

.656 

.660 

.739 

.740 

.803 

.803 

.8 

.816 

.817 

.843 

.844 

.870 

.870 



Eigenvalues for - 



A = 

1.4 



A = 

2.0 


a R 

Cl 

Co 

Cl 

c 0 


Thick-slab 

Present 

Thick-slab 

Present 

Thick -slab 

Present 

Thick-slab 

Present 


approximation 

solution 

approximation 

solution 

approximation 

solution 

approximation 

solution 


(eq. (J19)) 

(eq. (J6)) 

(eq. (J19)) 

(eq. (J6)) 

(eq. (J19)) 

(eq. (J6)) 

(eq. (J 19)) 

(eq. (J6)) 

0 

0.142 

0.132 

0.808 

0.808 

0.508 

0.508 

0.885 

0.885 

.2 

.252 

.247 

.820 

.820 

.540 

.539 

.891 

.891 

.4 

.403 

.402 

.836 

.836 

.585 

.584 

.898 

.898 

.6 

(*) 

(*) 

.860 

.860 

.656 

.656 

.909 

.909 

.8 

* 

(*) 

.900 

.900 

(*) 

(*) 

.930 

.930 


* a R > *>1 ( or ^1 * n branc ^ cu ^ I'oi)* 


with infinite reflectors have been recently computed by Kowalska (ref. 24) for a number 
of combinations of Cg^jj and c re ft ec t or . Some present results Co for °R ¥ 0 can 
be compared with the critical slab half -thicknesses of reference 24. The parameters 
are given in terms of C 0 and the present input quantities a D , and A by equa- 

tions (122). Figure 6 gives a few present cases (circles) for which Cg]^ was close to 
some of the points (diamonds) of reference 24; no attempt was made to compute exactly 
these points. The present cases for c s j a k ~ 1.11 are from A = 2 and 1.4 in table III. 
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Mitsis 


Critical slab half-thickness (ref. 2) I 

■ Critical sphere quarter-diameter (ref. 22)] 

O Critical slab half-thickness _ 

Present and Bowden 

$ Critical sphere quarter-diameter > (ref. 12) 

P O and 9 to scale of figure 


B 


□ 



Figure Critical dimension of bare systems. 



12 


□ Mitsis (refs. 2 and 22) 
O Kowalska (ref. 24) 

O Present 
Superscript is value of 


A 1.007 


iH 

CO 


•H 

■p 


aJ 

fH 


.01 

O 

1.01 


6 1.026 


1.01 

O 


1.01 

O 


1 . 01 ? 

o 


1.01 

o 

Oi.QX 


l 1 - 10 1.11 1.11 

6 ^ o O o 

1.13 1.12 1.11 


1.10 

O 


o 

1.11 


1 


1.11 

o 

l.n 

O 


0 .2 .h .6 .8 1.0 


c reflector 


Figure 6.- Critical half- thickness for finite slabs with infinite reflectors. 


Sample Reflected Slab Results 

The remainder of the results have been computed for A = 5. For a bare slab with 
A = 5, it can be seen from table I that there are five eigenvalues. The behavior of these 
eigenvalues has been studied as departs from zero for several values of Oq. In 

figure 7, results are given for a D = 0. The calculations show that the largest eigen- 
value £q is present up to a R = 0.9999. Apparently, this eigenvalue remains up to 
Or = 1, which is only obtained for C 2 < 1. All other eigenvalues disappear into the 
branch cut of Vqi at £ n = a R , labeled with an asterisk, which corresponds to a time- 
decaying mode, regardless of the value of C 2 - An asterisk is used in figures 7 and 
9 to 11 to indicate the points at which an eigenvalue or pseudo -eigenvalue coincides with 
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1.0 



Figure 7*- Dependence of eigenvalues on cr R . a D = 0; A = 5- 


VW\A Branch cut of v 
O Eigenvalue 

Q Pseudo-eigenvalue 





C 


Figure 8.- Relative location on the real £-axis of the branch cuts of Vq^ 
the eigenvalues £ , and the pseudo-eigenvalues for cr-^ =0.2 as read 
from figures 7 and 9 to 11. 
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the branch points of z^p. Even though such points appear to have a discrete eigenvalue 
type of time dependence, it is felt that they are properly part of the branch-cut integral 
contribution. Note that the branch points of are located at £ = -op> and 

£ = -Op) + Oj^ and that the limiting form of the condition which determines whether such 
points are eigenvalues (or pseudo-eigenvalues) no longer depends explicitly on o R or 
op). (See appendixes J and K.) The theorems of Lehner (ref. 8) apply for o R = 0 in 
figure 7. 

In figure 8, the manner in which the ^Q m branch cuts overlap and the location of 
the eigenvalues and pseudo-eigenvalues for cr R = 0.2 in figures 7 and 9 to 11 is indicated. 

In figure 9, results are presented for op) = -0.65 + 0.5opp. These results typify 
results for -op) values in the range between zero and (£ 0 ) 0 - _q> where the notation 

(^n) a Q means bare-slab eigenvalue, which depends on C 2 , an d a - The open and 

Rf 

closed circles represent eigenvalues as in figure 7 whereas the half-closed circles are 
pseudo-eigenvalues corresponding to s < -o min = -cr-p . Again the largest eigenvalue 
appears to remain provided that C£ > 1. Here, as in the next two figures, results for 
Ojp = 0 agree with the theorems of Hintz (ref. 9) which apply only for cp = 0. Basically, 
his result is that the strip Re(£) between 0 and -ap> belongs to the continuous spec- 
trum and that the bare-slab eigenvalues lying in this interval are not eigenvalues of the 
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r 


slab surrounded by purely absorbing media. He finds that there are no eigenvalues if 
-cr D > (5 0 ) a -o’ c * oes n °t di scuss the physical significance. It is seen from equa- 
tion (119) that for such cases, < l/c 2 and corresponds therefore to a time-decaying 
mode. In other words, stationary (critical) or time-increasing modes cannot disappear 
into the continuous spectrum as material properties are varied. In fact, when Or 4 - 0, 
such modes could not disappear into the branch cut of either. In figure 10, results 

are given for -Op + °r = 1 which, it may be remembered, implies C2 < 1. For this 
case, all the bare-slab eigenvalues lie in the continuous spectrum of reference 9 when 
°R = 0* In b otl1 fig^es 9 and 10, s = _cr m i n corresponds to £ = -Op. Figure 11 shows 
the behavior of the eigenvalues for ap = 1 and it is similar to that of figure 7. For 
°R = continuous spectrum of reference 9 lies in the strip -Op = -1 < Re(£) < 0. 


Here s = -a, 


mm 


corresponds to £ = 0. 


Eigenvalue and branch- 
point same 

Bare-slab eigenvalue 



Figure 10.- Dependence of eigenvalues £ n on cr-^. cr^ = cr^ - 1; A = 5 . 
Note that Cr> < 1 for this figure. 
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Summary of Numerical Results 

All numerical results indicate that real time eigenvalues for material 

reflectors are finite in number and tend to eigenvalues previously obtained for a vacuum 
as cr R — 0, as do the pseudo -eigenvalues for s < If the set <^ n ^> is empty, the 

neutron density is necessarily decaying in time. Conversely, if the neutron density is 
stationary or increasing in time, the set {?n} n °t em Pty- One also expects that 
when C£ > 1, a critical thickness should exist; that is, the largest eigenvalue £q must 
be present for large enough slab thicknesses for the given C£. This effect can be seen 
from table I. For example, if -or = 0.8, the eigenvalue for A = 1 is not present, 
whereas that for A = 5 would be, and represents a mode whose amplitude increases 
exponentially with time for C 2 > 1/0.975. That for A = 20 needs only C£ > 1/0.998 in 
order to represent a time-increasing mode. 


As pointed out at the beginning of this section, some speculations concerning the 
behavior of for reflected slabs as a function of the slab half -thickness a can be 

made, that is, if c m and o m are known what can be said about as a function 

of a. The following conclusions are based on the observation that if £q at Or = 0 
lies to the right of -Or, it appears to remain to the right of -or + or as or 
increases until -Or + Or = 1. (See figs. 7 and 9.) The dependence of on 

slab half-thickness is given in tables I to III and many more points are given m refer- 
ence 12. First, if -Or + Or = 1 (this inequality implies C 2 < l), the set is 

empty for all a. However, there may be pseudo-eigenvalues if -Or > 0. Next, if 
-o D + Or < 1, then two cases arise, depending on the value of or 


(1) When -Or > 0, then regardless of the value of C 2 , one can find 
an a* such that a < a* implies that the set is empty, whereas 

a > a* implies that the set is not empty. The number a* is 

obtained from the bare-slab result (Co) CT _q as 


|?o(«2.°a.a11 %=0 = -c D 


(123) 


(2) When -Or 2 o, the set is never empty. Thus, given c m , 

o m , a, and the bare-slab eigenvalues corresponding to C 2 , ovj, and a, one 
can say whether the set is empty. Furthermore, the number of eigen- 
values (Kx\) will not exceed the number of bare-slab eigenvalues _gj 

which are greater than -or. Finally, the number of real reflected-slab 
eigenvalues and pseudo-eigenvalues does not exceed the number of bare-slab 
eigenvalues. 
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COMPARISON WITH OTHER ANALYTICAL SOLUTIONS 


General Form 

It has been shown by using Case's method that the solution of the initial-value prob- 
lem of monoenergetic neutrons migrating in a finite slab (properties 03 , with infinite 
reflectors (properties can be written in the form 


= ^ U (x,/J.,t) + > 


s=s n 


p-^i-ci) 

+ 

2tu . 

)_a . 
u mm 

+ _L 

p- CT min+ ico 

27Ti . 

_0 min -ico 


e s ^ ds 


eSt ds (" a min < " c l) < s n) ( 124 ) 


In this equation, t is the real time multiplied by the constant neutron speed, <? m i n is 
the minimum of U]_ and 02, and each xp function is the sum of its definite parity 
parts x{/ ± . Some terms of the solution (124) will not be present if -o , m in^" cr l( : *- - c l) 
That is, if -aj(l - cjJ < -<7 m i n , then the branch-cut integral does not appear. Likewise, 
if all s n < -oi(l -ci), there are no residue terms. These discrete eigenvalue terms are 
characteristic of a finite slab whereas the branch-cut integral term is typical of a semi- 
infinite medium. The term ’Ji r u (x,/i,t) describes the behavior of neutrons from the ini- 
tial distribution f (x,/j.) which have not suffered a scattering collision and its definite 
parity parts are given in equations (104) to (107). The discrete eigenvalue terms in 
equation (124) are given by equation (117) whereas the integrand of the branch-cut inte- 
gral is given by equations (115) to (116). The definite parity parts of the last integrand 
are given by equation (103) and equations (70) to (73). The eigenvalues can be 

computed as was demonstrated in the previous section; thus, all terms in equation (124) 
can be calculated. 


Special Cases 

In all special cases of the present problem which have been solved using the Lehner- 
Wing technique (refs. 6 to 9), c^ = 0. In these cases, there is no branch cut due to 
Vqi(s); therefore, the branch-cut integral is not present in equation (124). It was shown 
that as ci — 0, the eigenvalues ^s n j> which are greater than -a m ^ n approach those for 
a bare slab as do the pseudo-eigenvalues for s < -a m ^ n . The solution if/ ± has the 
proper behavior as c^ — 0 since those terms of equations (65) and (66) which appear to 
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blow up in such a limit actually cancel when the contour C' is collapsed on to the por- 
tion of the branch cut of f2 m (z',s), 0 = z' =1. When the uncollided term is combined 
with the last integral, it is then seen that the solution (124) and the eigenvalues <{s n } 
have the behavior required by the theorems of Lehner (ref. 8) and Hintz (ref. 9). The 
present problem reduces to those considered by Lehner and Hintz, respectively, when 



and 





a l * a 2 ] 


j 


(125a) 


(125b) 


In order to describe the same physical problem in the slab as that solved by Lehner 
and Wing (refs. 6 and 7), one must not only have 


cj = 0 
oj = 0 


(126) 


but also 


f(x,ju) =0 (x < -a, ii > 0 and x > a, p. < 0) (127) 

In other words, neutrons from the initial distribution outside the slab cannot impinge on 
the slab faces at times t > 0. Lehner and Wing solved the time -dependent problem with 
boundary conditions: 

*(±a,M,t) =0 (jut % 0, t > 0) (128) 

Restrictions (126) and (127) in the present solution make l 2 ± (|U,s) and therefore 
A 2 ± (p.,s) depend only on slab properties. Then, in looking for solutions inside the slab 
(|x| < a), the inversion contour along Re(s) = -cr m i n can be deformed back to 
Re(s) = - 02 , and one picks up a residue contribution from any pseudo-eigenvalue in the 
region and thus obtains the Lehner-Wing results. That is, the solution has the proper 
form and all bare-slab eigenvalues are recovered. 
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The analogous problem for 
to the left of Re(s) = -CT m ^ n for 
that is, if 


Cj 4 0 in which the inversion contour can be deformed 
|x| < a is obtained when > cr-^ and f^(x,/x) s 0; 


f(x,p) = 0 ^|x| > a and 02 > (129) 

then all terms in I 2 ± (m,s) which contain s + factors in the exponentials are identi- 
cally zero and this allows the contour along Re(s) = -cr m i n to be deformed back to 
Re(s) = -02 when |x| < a. Such a deformation is not possible for |x| > a; for this range 
of x, one must stop at Re(s) = = -o^. If there are pseudo-eigenvalues in 

-02 < Re(s) < -0| = -a min (see, for example, fig. 9), they will appear in the solution for 
|x| < a as residue terms which have the exponential time dependence. They are however 
not eigenvalues for the reflected slab since such terms do not appear for (x| > a. 
Erdmann (refs. 15 and 16) solved the time -dependent problem for two semi-infinite media 
where an isotropic pulse of neutrons was introduced at the interface, and found that the 
inversion contour for x e medium m could be deformed to the left as far as 
Re(s) = -a m . In the present problem, such deformations can be made only when condi- 
tions (129) are satisfied. It appears that the contour Re(s) = -Omin cannot be deformed 
to the left of Re(s) = - 02 , since the implicit equation which determines A 2 ± (/x,s) (see 
eq. (13)) requires Re(s) = - 02 . Apparently, Re(s) = -02 is the edge of a continuous 
spectrum in all cases for the reflected slab. 

CONCLUDING REMARKS 

The present solution has been shown to have the required properties in all special 
cases which have been solved previously by others using the Lehner-Wing technique. 
However, in all these rigorous solutions, there was no scattering outside the slab. It 
was seen that with infinite reflectors on the slab and neutrons anywhere outside the slab 
initially, it is possible for some neutrons which have spent their entire history in the 
reflector to impinge on the slab faces at later times. Such neutrons have a collision rate 
which is characteristic of reflector properties and this condition, in general, places a 
restriction on how the inversion contour can be deformed in the transform plane. Two 
cases have been illustrated in which a further deformation is possible for the solution 
inside the slab by eliminating neutrons outside the slab initially, which can later impinge 
on the slab faces. This condition is equivalent to a further restriction on the Hilbert 
space which has been used in some of the above-mentioned rigorous solutions. The exact 
eigenvalue condition has been obtained, and real time eigenvalues have been calculated 
for a number of combinations of material parameters. The largest eigenvalues have 
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been shown to agree with the criticality results of others. Calculations also show that 
eigenvalues can disappear into the branch cut or continuum as material properties are 
varied and it was pointed out that all such disappearing eigenvalues correspond to expo- 
nentially time-decaying modes since the number of secondary neutrons per collision in 
the reflector was taken to be less than unity. It is expected (but has not been shown) that 
there is no drastic change in the shape of the solution when this situation occurs; one of 
the integrals in it probably has resonance-like terms caused perhaps by zeros of the 
eigenvalue condition on the next Riemann sheet. The assumption that the eigenvalues are 
real for arbitrary slab half -thicknesses has been made. This assumption has been shown 
to be true for thick slabs in this report and it has been proved rigorously by others for the 
above-mentioned special cases. On the basis of sample calculations, it is concluded that 
if one is given the material properties as well as the bare -slab eigenvalues corresponding 
to the slab properties, then he can conclude whether there are eigenvalues and the maxi- 
mum number of them. 

Perhaps the present results can serve as a guide for a rigorous Lehner-Wing type 
analysis of the reflected-slab problem. If the eigenvalues are all real, one might be able 
to prove it in such an analysis of the present problem. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 12, 1970. 
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APPENDIX A 


SUMMARY OF ELEMENTARY SOLUTION PROPERTIES 

In this appendix, the elementary solution properties derived by others (refs. 12 to 
to 16), following the lead of Case (refs. 2 and 11), are summarized. These solutions are 
obtained from equations (21) to (23) and are given by equations (24) to (26). Such solutions 
are complete and orthogonal in the following sense. A function, say g(ju), satisfying very 
weak restrictions (see, for example, appendix G of ref. 2) for -l=a! = p.S/3 = l can be 
expanded as follows: 

(1) Full range (a= -1; /3 = 1): 


g(/i) = [a m ^ 0m (p) + b m^ 0m (M) 


(s) + f A m {v) s) dn (Al) 

J -1 


where the notation 6 m (s) was defined by equation (30). The orthogonality relations used 
to determine the expansion coefficients in equation (Al) are 

(MjS) dp. ^ ^ A nl (p) s) dn — A m (n ) v {v ,s) {v ,s) 


and for s e S mi 


J s) (p (m) dn = 0 

f W Vn (M) (M) d/J. = 0 

J-1 Tm 0m 

y 1 ^ = 2 Cm<Tffl!; 0m ^m( ±ll 0 m ,s ) 


(A2) 


where 


n m(W s ) = i n m^ s ) 


(A3a) 


z=n, 


0m 
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and 


+ . . . iffc rn a rn ^ 

(M = A m (y,s) ± 


(2) Half range (a= 0; /3 = 1): 


g(fj.) — a. m o7^ (fj.) 5 m (s) + C dv 

0m 


Here the orthogonality relations for s e S m j are 


pi pl 

J 0 W m M J 0 A m( u ) d,/ = A m (i/') W m (^) fi m V,s) a^(v',s) 

w m (^) ^rm.(tt ; s) ^^(pt) dpt = 0 

w n>> %>' s > ^-. 0m (ft) = t ’ c m CT m ^0m X m ( _i/ 0m> s ) 

J q W m ( M ) <?t, Qm (M) ^±v 0m (M) dpt = T ( C - g f' - m ) 2 X m (±y 0m ,s) 

w n» ^m(-^)(h,s) dpt = ^ Qm X m (~v,s) 

c (j 

J 0 W n>> ‘Pmv'M <P m (_ v )(M,s) d M = + ^’)X m (-v,s) ^(.^Ks) 

pl 1 

J W n>) dh = 2 c m a m^ 


where 

W m (M) = 


C m° r m/ J ' 


2fi m (~s) p 0m + 


(0 = M 


(A3b) 


(A4) 


(A5a) 


< i) (A6a) 
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and 

x rc(M,s) _ 

(MjS) 


(0 = M = 1) (A7a) 


with X m (z,s) given by 



(A8a) 


Several other equivalent expressions for X m , generally referred to as the X-identities, 
are given by 


r 


ft m (z,s) 


X m (z,s) =/ 


(^Ora ‘ z 2 )^m(°° 5 s ) x m(-z,s) 


c m° r m C 1 

^0 (v. 


M dja 


l v 0m ' V 2 ) X m(~V’ s ) (M - z) 


(A9a) 


The orthogonality relations for s e S me for the expansion (A4) are 

r l r l . 

\ W m (n) dji \ A m (v) di/ = A m (t' 1 ) W m (i/’) JZ m (v’,s) Sl m (v',s) 

^0 J 0 

nl i 

J 0 W n>> = 2 Vm^'W'bs) 

pi j 

W m (M) ( Pmv^ ,s ^ ^ = 2 c m°m 1 ' 
where 


) (A5b) 


W m (M) = 


c ra ff mM 


2^m(°°) s ) x 0m( 


(0 s M = 1) (A6b) 


and 
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XQ m (M-jS) ^ (MjS) 

x 0m^> s ) n m “(^ } s) 


with X Qm (z,s) given by 


^Omf z > s ^ 



and the identities 




n m (z,s) 


X 0m( z > s ) ~ \ 


s) X 0m (-z,s) 


+ c m G m (*' 


.in -ixi i d/i. 

2£2 m (°°,s) J 0 X 0m (-ju, s) (m - z) 


(0 § p s 1) (A7b) 


(A8b) 


(A9b) 


These half-range orthogonality relations and identities are obtained by extending the 
time-independent results of Kuscer, McCormick, and Summerfield (ref. 20). 

A result, due to Kuscer and Zweifel (ref. 14), which is needed to continue solutions 
analytically follows from equations (A9a) and (A9b). For a fixed value of z, X m (z,s) 
does not become XQ m (z,s) as s crosses C m . However, it can be shown from the 
first line in equations (A9a) and (A9b) that 


( v 0m " z)x m (z,s) 


s "* C m 
seS mi 


x 0m( z ’ s ) 


s— C 
seS 


m 

me 


(A10) 


By use of reference 14, the function XQ m (z,s) is redefined to be continuous as s — C m 
by equation (60). Such a function of the two complex variables z and s has the fol- 
lowing analytical properties (ref. 14): 

Fixed s: no singularity in z-plane cut along (0,1); one simple zero at z = z'om^’ 
Re M= 0) ° nlyif seS mi 

Fixed z: no singularity in the s-plane cut along (-cr m , -cr m (l - c m )j; one simple 
zero at s = -cr m + c m a m z tanh - -*- — for Re(z) > 0 

Note here that XQ m (z,s) is a nonvanishing analytic function of z and s for Re(z) < 0 
and s 4 (-a m , -a m (1 - c m )), the branch cut of i^Q m (s). 


60 



APPENDIX B 


DERIVATION OF ^ mp± (x,/i,s) 


In this appendix, explicit forms of ^ mp± (x,p,s) are obtained. Consider, for 
medium m, the function g m (x,/i;xo) as 

-D m (x„) e ( s+cf ra)( x-x o)/ y Om 6m(s) 

- C m (x 0 ,x) s) e'( S+<T m)( x - x O )/" dv (x < xq) 


g m (x,ii;x 0 ) =< 


(Bl) 


C m( x 0 ),„ 0m („e-( SWm )( !t - X °)/^6 m (s) 


(x > x 0 ) 


The expansion coefficients in equation (Bl) are to be determined so that g m ^x,p;xQ^ 
satisfies equations (42) and (43); that is, on putting the expansion (Bl) into equation (43), 
in the limit x — xq 


M 


= j^-'m( x o) <P^ m (M) + D m( x o) ‘^om^ + C m( x 0^) 


(B2) 


This equation is a full-range expansion (see eq. (Al)) of the function f m (xg,p,yM and 
use of the orthogonality relations (A2) gives the coefficients as 


C m ( x 0’ 1 ') - ‘ , f . f m(*0>») <U* 

v^(v,s) n m (v,s) J -1 v 


and, if s e S 


mi 


C m( x o) = f « f m( x 0^) <Pv n 0*) d M 

V 7 c m a m , 0m 2fi^ 0m ,s) d -l ' Om 


(B3a) 
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and 


Dm(Xo) = c m ^A(- W) t dtl 


(B3b) 


However, one needs expansion coefficients for f m± (xo,p.y/i- It follows from equation 
tions (B3) that 


C m±( x 0> 1 j s ||Cm( x 0^) + C m 


and, if s e S mi 


nu> 


and 


Cm±( x 0^0m) = ifm( x o) T D m (-xo)] 


C mi( x 0>' l '0m) - ■|[l ) m( x o) + c m(" x o)J 


(B4a) 


(B4b) 


are the expansion coefficients of that is, equations (47). 

In order to construct ^ m p±(x,p.,s) according to equation (41), note that for m = 2 


T (...) dx 0 = P (...) dx Q = T (...) dx Q + P (.. .) dx ( 

J medium 2 J -a J -a J x 


(B5) 


Upon using equation (Bl), one obtains i// 2 p(x,p.,s) as 
^ 2 p( x ,M,s) = ft* C 2 (x 0 ) e ( S+a2 ) X °/ I '02 dxc 


^02 (X ’ M,S ^ 


f D 2 (x 0 )e-( S+ ° 2 W 1 '02 cb , 


i// 2 i;(x,M,s) dv 


't / _ VQ2 &,V; s ) S 2 (s) 


^2(-y)( x ’^» s ) ( B6 ) 
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The definite parity particular solution \J/ 2 ^ ± (x,fJ.,s) ^ en obtained by using equa- 
tions (B6) and (B4) as 


^2p ±( x >M,s) = 


if C 2± (* 0 ,*0 2 ) e^^ 02 

£ 0 2± ( Xo ^ 02 ) 




*_ VQ2 b,ll,B)) 6 2 ( s ) 


^2 j/( x ^^ s ) ^ 


i! r c 2 ±( x o> ! ') e (s+ ^ x ° /i ' dxo 


^2(-!/) (x ’^ s) 


(B7) 


That equation (B7) is a solution of equation (18) for m = 2 can be seen by direct substi- 
tution as follows. The \p mv (x,ii,s) in equation (B7) are solutions of equation (21), the 
homogeneous equation corresponding to equation (18). However, their coefficients in 

equation (B7) are functions of x so that some additional terms are obtained from the 
3 

— operation. Thus, 


M|jC2 ± ( x ^ 02 ) + c 2±( x } -^02) <P- 

+ J i C 2± (x,y) cp 2v {^,s) dJj = f 2± (x,/i) 



(B8) 


which is an identity since according to equation (47), the C 2± are the full-range expan- 
sion coefficients of f 2± (x,/i.)//i. 

To get if/ip(x,n,s) according to equation (41), first note that 


r (...) dx n = r (...) dxg + r (...) dx n 

^medium 1 U J -~ J a u 


7>x p-a p 00 

\ (...) dx 0 + \ (...) dx Q + \ (. . .) 

J-OO 


dx n (x < -a) 


-< 


(B9) 


p-a px P°° 

\ (...) dx 0 + \ (...) dx 0 + \ (...) dx 0 (x > a) 
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By following the same procedure as before, one gets l / / ip ± ( x >M> s ) as 


^l p± (x,p,s) 


J ^ C l±( x 0^0l) e( S+Cri ) X °/ Z ' 01 dx Q ^ q1 (x,M,s) 6 1 (s) 


4/ lv {x,n,s) dv 


-J x c l±( x 0.- 1, 0l) e ( S+ffl )V 1 ' 01 dx. 


^_„ oi (x,/i,s) 5 1 (s) 


J c l±( x 0^) e( S+ai ) X ° /" dx c 


di/ (x < -a) 


and 

^lp± (x ’^’ s ) = 


£ C 1± ( V01 ) e^Aoi ^ 

nl n-a. , ('s+cr 1 '\x n /y 

^2 c u*o--")°' {s+ ° l)xo/v «* 

£ c l±^.*t)x) e( S+<7l ) X °/ ! ' dxo 


^ 01 ( x >M,s) 6 1 (s) 


4/ lv (x,iJ.,s) dv 


)_oo C 1 ±( X 0 ,I/ ) e 


( s+CT i) x o/ y 


dxr 


>/ / _ I , 01 ( x > M,s) 6 1 (s) 


^l(-j,)( x ,M,s) dv (x > a) 


(BlOa) 


(BlOb) 
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Again, it is easily shown by direct substitution that equations (BlOa) and (BlOb) are solu- 
tions of equation (18) for m = 1. By introducing the F functions of equation (46) and 
allowing x to take on negative and positive values, it follows that equations (B7), (BlOa), 
and (BlOb) can be written as equations (44) and (45). 

Also note here that the C m± coefficients of equations (B4) have the property 


C m±(- x 0>-^) = : F C m ±( x 0> 1 ') 

C m±(- x 0»^0m) = + C m±( x 0^0m)J 

so that it then follows from equations (46) that 


(Bll) 


F 2± (a,-w,s) = +F 2± (a,w,s) 
F ± (-a,-w,s) = TF ± (-a,w,s) 


(B12) 
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APPENDIX C 


TWO-MEDIA FULL-RANGE EXPANSIONS AND 
ORTHOGONALITY RELATIONS 

In this appendix, some results of references 15, 16, and 20 are summarized and 
extended. In reference 15, it is shown that a function, for example, h(p), satisfying very 
weak restrictions for p on the interval -1 ^ p = 1 can be expanded as 

pi p0 

h(p) = ajoJj , (p) S-^s) + b 2 cp_ v (p) 6 2 (s) + \ A^y) ^ (p, s) d^ + \ A 2 (v) <p 2l ,(p,s) ^ 

U1 0 Z Jq ^-1 


This equation is a two-media full-range expansion of the function h(p) and the expansion 
coefficients in it can be determined by using orthogonality relations which are easily 
determined from the time-independent ones of reference 20. For 6^(s) = 6 2 (s) = 1, that 
is, s e S^ n S 2 p these relations are 

r l pl "I 

\ W(p) $j,,(p, s) dp \ A(v) $j,(p,s) dv = A(v') W(p') Q, + (v',s) 12 {v',s) 

J_1 J - i 


W(p) $„(p,s) <Pi/ 01 (p) dp = 0 
J W(p) $ y (p, s) <p_ j/ Q2 (p) dp = 0 


f W(p) (p) cp_ (n) dp = 0 
J-l V 01 02 

pl 

J W(p) <P_j, 01 (p) dp = pc(v) cr(t’) ^ 01 ^ 0 2 " "OlJX^Ol* 8 ) 

pl 

J $ p(M,s) ^(m) dp = vc(v) o{v) v Q2 (p 01 - ^X^OB* 8 ) 

W ^) ^ Q1 (fi) dM = -(~ T 01 ) fal * v 02)*(±V01>s) 


Equations continued on next page 
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^ W ^) (p ±VQ 2 {li) ^ = C 2 T° 2 ) ("02 T v oM ±V 02’ s ) 
I_! W(M) %1 (m) % 2 (m) dl± = ~2 C ^ C ^ V 01 P oh{ V 02’ s ) 

pi -i 

J W ^) ^-^ 02 (m) dM = 2 Cl ° r l C 2 a 2 I/ 01 W(- y 0D s ) 


c(v),o(v) = 


* v (v-,b) = 


c l’ a l 
= J 

(v > 0) 

[ c 2> a 2 

{v < 0) 

|S>i„o*,s) 

O' > 0) 

\ 

<P 2 >> s ) 

V 

( v < o) 

(flffos) 

{v > 0) 

1^2 (v,s) 

(KO) 


QT(v,s) = 


X(z,s) = X 1 (z,s) X 2 (-z,s) 


l p 02 + 0 X 2("^ S ) W l^) 

-(^01 " ^) x l(^s) W 2 (-v) 


( v > 0) 
( y < o) 



All remaining quantities have been defined in appendix A. 

Rather than write out explicit orthogonality relations for the other three regions of 
the transform plane in the present notation, one introduces a function which is continuous 
as s C m . From the results of reference 14 quoted in appendix A, it is seen that one 
such function is given by equation (59) and can be written by using equation (60) as 
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7^02 " Z ) X 2( Z > S ) 

(s e s n n s 2i ) 

(i/ 0 i - z)X 1 (z,s) 

(^02 " z ) x 2( z »s) 

( s 6 s le n S 2i) 

X 0 i( z ,s) 

x 02( z > s ) 

(s e s li n s 2e ) 

(^01 ' z ) x i( z > s ) 

x 02( z ? s ) 

x 0 i(z,s) 

( s 6 s le n S 2e) 


In terms of this function, W(^) can be written as 


(C4) 


(°°s) 


X 0 (-y,s) 


W(^) =< 

C2^2 V 1 

^2£2 2 («,s) Xq(v,s) 


(v>0) 


(C5) 

(v< 0 ) 


The function x( z > s ) is expressed as 


X 0 (-z,s) 

%(z,s) 

(^02 + z ) 

(i/qi - z)%(°°,s) 


1 Q2( g » s ) 

(v 0 i - z)X 0 (z,s) (V 0 2 + z )^2(“ s ) 


(Re(z) > 0) 


(C6) 

(Re(z) < 0) 


In order to obtain a two -media expansion in the form of equation (Cl), one generally has 
to switch some continuum solutions in one medium to those in the other. From the 
explicit form of cp mv (ix,s) (eqs. (24)), it follows that 


c l a 1^2^> s ) " c 2 a 2^1^> s ) =k6(/i - v) 


where 


k = 5(0^ - c 2 a 2 ) + - c 2 ) 
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It is seen that when the two media are the same, k = 0. This quantity can be expressed 
in a number of different ways, and several that will be used are 


'c 1 a 1 X 2 (i',s) - c 2 a 2 X 1 (^,s) 
k =/c 1 a 1 ft 2 (y 01 ,s) 


(C8) 


j c 2 a 2^l(Vo2> s ) 


The orthogonality relations (eqs. (C2)) can now be written in terms of Xq(z,s) and 


k as 


J 1 W(ju) ill, s) dp, ^ A{v) s) dv = A(v') W(v') n + (v',s) Q~(v', s) 

pi 

j W(ju) * v {n,s) <p VQ1 (li) dji =0 

pi 

^ W(m) dn - o 

W(n) d " ’ 0 

P yc(l^) cr(v) iM-jk 

li W(fl) *>’ S> d(i ■ 4tt 2 Ks)X 0 (-*01,sUv 01 *-) 

, fe-m-ci)’ SJgjVoi'*) 

[,»W V„ 01 <M) <»(*=< 


} (C9) 


c l g 1^01 k 

8^ 2 («>,s) Xo(-i/ 01 ,s) 


Equations continued on next page 
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<P. 


( m ) <p± Vn M d ^={ 


1> Q2 ±V Q2 


c 2 a 2^Q2^o(-^Q2> s ) 

8^1 (°°, s) 

[1 >,2 ^2 (^ 02 ’ s ) 
'{| C 2 a 2^02) 


x 0(-^02> s ) 


<P Vq1 M = 


c l g 1^01^02 kX 0(- t/ 02>s) 
4%(~,s) (v 01 - ^ 02 ) 


^_jy 02 (M) <P-v 01 (m) d/i = 


c 2 cr 2 I/ 01^02 k 1 

4fi2(°°> s ) (^02 " ^Olj X 0( -J/ 01 >s ) 


J 


These expressions appear to be more complicated than the corresponding ones in equa- 
tions (C2); however, the orthogonality relations needed for all regions of the transform 
plane are given by equations (C9). That is, for s e Sj e n S 2 i, the proper orthogonality 
relations are the first, third, sixth, and eighth equations of equations (C9) with Xq(z,s) 
given by equation (C4). Note here that Xq(z,s) always appears in equations (C9) with 
Re(z) < 0. It follows then from appendix A that for Re(z) < 0, Xg(z,s) is a nonvanishing 
analytic function of both z and s except for the branch cuts in the s -plane due to 
i/ 01 (s) and v Q2 (s). 
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THE TWO-MEDIA FULL-RANGE EXPANSION FOR THIS PROBLEM 


In this appendix, it is shown that application of the continuity condition (eq. (20)) 
results in a two-media full-range expansion of the type discussed in appendix C. For 
x = a, one readily obtains from equations (31) and (20) upon using the explicit forms of 
^mc± anc * ^mn+ S iven by equations (32), (39), (44), and (45) that 


mp± 


0 = a, 


|^ 02 (a,/X ’ S) * ^ 02 (a ^ ,s) 


£ A 2* W 


5 2 (s) 

^ 2 v {a.,ji,s) ± ^ 2 (-i/)(a,M,s) 

a l±^V 01 ( a >^> s ) M s ) + A l±(- p ) ^liy( a »M,s) dz^ 

pi 

+ F 2± (a,^ 02 ,s) ^i/ 02 (a»M,s) 5 2 (s) + ^ F 2± (a,i/,s) ^ 2y (a,M,s) d* 

- F 1± (-a,v 01 ,s) ^ 01 (a,M,s) ± ^-i; 01 (a,M,s)| « 1 (s) 

- J F 1± (-a, v,s) |^iy(a,jLi,s) ± ^i(_ y )(a,M,s)] dz/ 


(Dl) 


It was indicated in appendix C that according to Erdmann (ref. 15), the functions cp.. in), 

‘ V 01 

cp_ v ( n ), <Pi v (ijl, s) (0 ^ v ^ 1) and <p 2 (ju,s) (-1 ^ y s o) form a complete orthogonal 
0 2 

set of basis functions for the expansion of h(/iz) (-1 S ju ^ 1) for s e S^ n S 2i . (See 
eq. (Cl).) However, equation (Dl) also contains terms in which <p 2l/ (/i,s) (0 ^ ^ 1) 

and <pij,(/i,s) (-1 ^ v S o) appear. These continuum solutions must be replaced by 
corresponding continuum solutions for the other media. One uses the relationship (C7) 
to do this; that is, 




<Plj;(ja,s) = 


— — cp 1 (jU.,s) + k 6(v - ii) 
c l a l Iv c l°l 


Cl O’ '| b- 

<P2v&’ a ) 


c 2 a 2 


H(z^) 


H(-z/) 


(D2) 
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where 


(v > 0) 


( v < 0 ) 


When explicit forms of the elementary solutions and equations (D2) and (D3) 
are used in equation (Dl), one obtains the two-media full-range expansion 


M^) = F 1± (-a 3 !/ 01 ,s) ± a 1± e ( s+<Tl ) a Aoi <p v fa) 6i(s) 


t e 


«,W 


+ < |Fi ± (-a,^s) ± A 1± (-r/)| e 


'( s+cr l) : 


F 2± (a,y,s) + A 2 ± (i/) e 


' ( s+0 2) a / 


(p lv fa,s) dv 


£ [2| Fl± (. a ,-,, s) a-( S+ °l) a /^ . a 2± (-,) (D4a) 


where hfa) is given by 


h(/x) = ^ H(m) §2±( a > M,s) + A 2± (/i| e 


•(s+CT 2 )a//i 


Fi±(-a,-M,s) e ( S+<T l) a / M 


[ F 2 ±( a ^ 02 > s ) + a 2 ±] e ( S+a2 ) a / V ° 2 ^( fx ) 5 2 ( s ) 


. . fs+CT-Aa/v n1 

T F l±(- a ^0D s ) ^ 01 6 1< S ) 


(D4b) 
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The orthogonality relations (eqs. (C9)) can be used on the expansion (D4) to obtain equa- 
tions which determine the remaining unknown coefficients implicitly. However, it is 
convenient to introduce first the E m± coefficients given by equations (50). One then has, 
after some algebraic manipulation, equations (51) to (48). 
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APPENDIX E 


CONTINUATION TO THE COMPLEX PLANE 


In this appendix, the way in which equations (50) to (58) are extended to the complex 
plane n — z with s considered a parameter is outlined. As shown in reference 21, 
functions such as those introduced in equations (50) are extendable. The particular 
grouping of terms in equations (51) to (58) indicates some integrals and residues which 
go together. 

The first functions to be considered are the F m± (x,w,s) functions given by equa- 
tions (46). In equations (55) to (58), these functions appear with Re(w) >0 so one con- 
siders the functions L m± (x,^,s) given by equation (69). When the explicit expressions 
of equations (46) and (47) are used, one can see that for f m± (xQ , m) extendable /i — z 
without singularities in the finite z -plane, then L m± (x,t',s) can be extended to 
L m± (x,z,s) given by equations (67) and (68). As z — v e (0,1), it can be seen that the 
limiting values of L m± , namely L m ±(x,i',s) and L m; ^(x,y,s), are identical. Thus, 
L m ± does not contain the branch cut of £2 m (z,s) as one might be led to expect from 
equation (67). There appear to be no other singularities of L m± in the finite z-plane, 
Re(z) > 0 and Re(s) > a m . It follows from equation (67) that 


V±( x, ^0m’ s ) " 2 c m a m v 0m^ m ( I '0m) s ) ■^m±( x > I ' om > s ) e ^ ^ ^ (s e S m ^ (El) 


In order to extend the functions I m± (f) to the complex z-plane, one needs the 
identity 

c l°'1^2~( i/ > s ) _ c 2°2 k 

Cgagft^jS) fi-f^s) c l a l c l°l c 2°2 


cia 1 X 2 (^,s) + c 2 tf 2 A 1 (i',s) 


«j>,s) (v,s) 


(E2) 


which can be verified directly. On using this identity, one finds that I m± (^), given by 
equations (55) and (57), can be written, respectively, as equations (65) and (66). The 
restriction Re(s) > on these equations comes from the fact that L m± for both 

m = 1 and m = 2 occur in each I m± . More will be said about this restriction later. 
The contours C' are given in figure 3. By letting z = Vq 2 in equation (65) and 
z = in equation (66), it can be seen that 
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l 2±{ v 02> s ) mJ 2±{ v 02) 


(E3) 


Thus, the inhomogeneous terms of equations (51) to (54) are seen to be extendable and 
related as shown in equations (E3). For z — z^ in l2± (eq. (65)) and z — z^ in 
Il± (eq. (66)), these functions might seem to be singular. However, upon examining the 
residues it is seen that this is not the case. Thus, the I m± (z,s) appear to be analytic 
in the finite z- plane, Re(z) > 0 and Re(s) > 

In equation (51), one now lets v — z and, for Re(s) > -a m ^ n and Re(z) >0 in 
the finite z -plane, finds that E2 ± (z,s) is given by the inhomogeneous term l2 ± (z), a 
term involving a2 ± if s e S2i and an integral over E2 ± (m,s) (0 = M = 1). A singu- 
larity occurs in the integrand when either &2"(/U,s) or (m 5 s) vanishes and this 
condition occurs for s e C2. However, for this case, it is seen that one obtains from 
equations (51) and (52) that E2 ± (i / q 2> s ) is related to a2 ± . It appears that E2 ± (z,s) is 
analytic in the finite z -plane, Re(z) > 0, Re(s) > -a min and can be written as equa- 
tion (63). By following the same procedure with equation (52), one obtains an equation 
which is easily seen to be equation (63) evaluated at z = v Qg,’ that is, E2 ± (z'q2j s ) anc * 
a2 ± are related as 


E 2±{ v 02> s ) “ | c 2 CT 2 t '02 fi 2( i '02’ s ) 


l 2± 


X s+o 2>h 


02 


(seS 2i ) (E4) 


In a similar manner, from equations (53) and (54) with v — z and from using equa- 
tion (E4), one obtains equation (64) and again it follows that Ei ± ^oi,s) and ai ± are 
related as 

E 1±( I/ 01» S ) = ^ c l cr l l n l (^0 1 » s ) a l± e ^ ^ / 01 ( s 6 s li) ( E5 ) 


It also appears that E^ ± (z,s) is analytic in the finite z -plane, Re(z) > 0 and 
Re(s) > -a min* 

The solutions z// mc± and i// m p ± can now be written in terms of the E m± as 
shown in equations (70) to (73). 
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INVESTIGATION OF THE ASSOCIATED EIGENVALUE PROBLEM 


In this appendix, the associated eigenvalue problem, that is, the problem for which 
f(x,/x) = 0, is considered. The inhomogeneous terms I m± (given by eqs. (65) to (67)) 
can be seen to be identically zero when f(x,p.) is zero everywhere. Solutions for 
I m ± = 0 are denoted with a bar, that is, i// m± . The unknown expansion coefficients 
for the associated eigenvalue problem E m± are given by equations (63) and (64) with 
I m ± = 0. It is seen from such equations that E m± can be determined only to within an 
arbitrary factor independent of z and that Ej ± depends on E 2 ± . Furthermore, the 
original normal -mode expansion coefficients for the eigenvalue problem are given by 
E m± (i±,s) (0 2 S 1, m = 1,2), E 1± (^ 01 ,s) (s e S n ), and ( s € s 2i)- 

Therefore, one must examine solutions of such equations as a function of the transform 
variable s for z — /a with the contour C' collapsed onto the branch cut (0,1) due 
to Q, 2 (z',s), and for z = VQ m when s e S m ^. This procedure is followed for all s in 
some right-half s-plane and it is convenient to divide the plane into three regions: 
s e S 2 e , s e S 2 i? and s e C 2 . 

When s e S 2 e , does not vanish within C' so that equation (63) with 

l 2 ± = 0 can be written as 

pi 

E2±(^,s) = ±J Q K (i±,u) E 2± (^,s) d^ (s e S 2e , 0 = /i = l) (FI) 

where 


K(iM 


k ^ 2 ^°°’ s ^ X 0 (-^ s ) v -2(s+02)a jv 

2 fii(°°,s) f2 2 +(i/, s ) fi 2 "(^s) (y + m) 6 


(OS^Sl, M = ^02) ( F2) 


When se S 21 , fi2( z '> s ) vanishes inside C' but not on the real interval (0,1). 
As C' is collapsed on to (0,1), a residue term appears so that equation (63) with 
I 2± = 0 takes the form 


E 2± (M,s) 


± 


k ^2(°°> s ) Xq(-m,s) X 0 (-^ 02 ,s) - 2 (s+a 2 )a/z^ 02 

c 2 a 2 %(°°> s ) ^ 2 ^ 02 ^) (^02 + M ) 6 2± 



± J K (n,v) E 2 ± {v,s) dv 


(s e S 2i , 0 < (1 < 1) (F3) 
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for z = ix where K (^,v) is given by equation (F2). However, equation (63) with l 2 ± = 0 
must also hold at z = ^g 2 anc * this requirement gives an additional constraint on solu- 
tions of equation (F3), namely, 


k ^K 8 ) X 0 (~ t/ Q2 ,s ) o - 2 ( s+CT 2) a A02 ^ 


2±l ° 2, ) " c 2*2 Ol(-,s) 2, 02 ^(, 025 s) 


2±( t ’02 ,s ) 


pi 

± J Q k (^ 02^) E 2±^’ s ) dv 


(s e S 2i ) (F4) 


When s e C 2 , the curve separating S 2 i and S 2e , ft^Xz^s) ^"(^s) vanishes for 
some v on the interval (0,1); that is, ^2 is real and lies on (0,1). By setting 
1^02 = V> one can P u t equation (63) with I 2± = 0 in the form 


E 2±(P> S ) 


1 k n 2 (°°,s)X 0 (-ii,s)X 0 {-Ti,s) " 2 ( s+a 2) a A ? , s 

± 2 c 2 a 2 fli(eo,s) 4 (t?,s) (t? + M) ® 2± ^’ S; 


k^2(°°> s K t ^ -n f 1 e 2 ±(*',s) X 0 (-^,s) -2(s+a 2 )ap v & v 

± 2 ^ 1 ( 0O ,s) 0 M,S y=r? J 0 fi+(i7,s) n 2 "(i/,s) 6 *' +fA 


(s e C 2 , 0 = n,n s l) (F5) 

Note that this equation would be obtained from either equation (FI) or equation (F3) for 
s — C 2 from s e S 2e or s e S 2i , respectively. 

For arbitrary complex values of s, the kernel K(/i,^) which appears in equa- 
tions (FI) and (F3) is not symmetric since 

K(fX,v) 4 [k(^,m.)J* (Im(s) 4 0) (F6) 


where the asterisk denotes complex conjugation. Note, however, that when Im(s) = 0, 
the unknown functions E 2 ± (ju.,s) can be redefined so that a symmetric kernel is obtained. 
Solutions of equations (FI) and (F3) depend on the behavior of K(i±,v) generally through 
the quantity b2(s) given by 


B2(S) = £ £ diti dv 


(F7) 
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To study B 2 , one introduces the nondimensional parameters £, Oq, and A given 

by equation (118) with £ = a + i/3. Note that a, p, cr R , ct-q, and A are real whereas 
and A are nonnegative. In terms of these quantities, one has 


K(p,i/) 


f[q(q R - 1) - a D ] 2 + /3 2 (a R - l) 2 ] - D 2 + P 2 ] |X 0 (- M;S )| 2 ^ 0 (-,,s)1 2 ,2 -4Aa/v 


4 j( a + a D " a R) 2 + £ 2 ] 


S?2(^,s) 

2 

n 2 ’(^s) 

c 2 a 2 


c 2 ff 2 


(v + (i ) 2 

(0 § p,v £ 1) 


(F8) 


To make estimates of the function |Xq(-P,s) |^, one uses an integral representation 
of the single-medium X -function given in reference 14, namely, 


X 0m(-^> s ) = ex Pte£ ico lo £e 


n m (z', s ) 

dz' 

3n(~ s )_ 

z' + pj 


(F9) 


Upon letting z' = iy and using S2 m (z’,s) = S2 m (-z’,s), one sees that equation (F9) 
becomes 


x 0m (-/ t > s ) = ex P<^ 


log. 


^mCly^s) 

dy_A 

nm(°°,s)_ 

y 2 + p 2 J 


(F10) 


which is real for s real. In terms of the quantities of equation (118), one finds 


2 

| X 0 i("M,s)| = exp 


L ft? +_?d -_ C R^ 

Y o e j (a + o D - c R ) 


& + a D ~ g Rg(y)| 2 + 0 2 dy 

) 2 + j3 2 y 2 + H 2 




(FI 1) 


= exp v J l0 & 


00 ' [a - g(yj] 2 + /3 2 


dy 


(a - l) 2 + /3 2 y 2 + p 2 


where 


g(y) = y tan -1 i 


(F12) 
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It follows from equation (F12) that 0 = g(y) = 1 for 0 ^ y ^ and that it is a monotonic 
increasing function of y on this interval. Furthermore, since 


2/xf dy _ 
17 ^0 y 2 + ju 2 


the following bounds (perhaps rather loose) are obtained for XQ m (-ju,s)|' 


ls i 

]j(a + a D ~ a R y + , 


(F13) 


(F14) 


for a > max(l, -ct d + OpJ, and 




I (<* - l) 2 


— 2 i |x 02 (-M,sf s 






I fl2 . l9 / / (a + a D ) 2 + /S 2 

,/7 — T2 2 = X 0l(“^’ s ) 2 = maxl >, 7 ; 

]/(a + a D - a R ) 2 + /3 2 \ ^ ( a + a D ~ a R. 


_) 2 + /3 2 




(FI 5) 


for a < max(l, -Op + and j3 ^ 0. Note that in the £ -plane, the points ? = 1 and 
£ = -Op + a R are the right ends of the branch cuts of 1^02 and v 01> respectively, and 
those cuts lie on the real £-axis. The left ends are at £ = 0 and £ = -Op, respectively. 


The functions 


^(^s) 


c 2°2 


are easily found to be 




Co CT< 


22 


= (a - v tanh - -*- vf + ^3 ± 


(FI 6) 


It was pointed out in the text that the curve C 2 (fig. 2) is given by 


a' - tanh ^ 

1 t 7 r 


(F17) 
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The parametric form of this equation is 


1*1 -? 


a' = v tanh“l v 


(0 ^ v $ 1) 
(0 ^ v z 1) 


(F18) 


It is seen from equation (FI 6) that 




c 2°2 


are the squares of the distance (in the 


£ -plane) from the point (a,/3) to the points (a'(v),+i 3'(u)), respectively, which lie on the 
curve C 2 . Since these functions appear in the denominator of |K(/i.,z')| the inte- 
gral (F7) will not be bounded when a and /3 are related as in equation (FI 7). One 
defines D m ^ n (o!,/3) as the minimum distance from the point (a,/3) to the curve C 2 for 
0 = ^ = 1; that is, 


~ min i 


Ci2(v,s) 


Co Oi 


2 U 2 


^2 


Co OV 


22 


(OSi/Sl) (F19) 


and D min (a,/3) ^ 0 for (a,/ 3) / C 2 . Therefore, from equations (F16) and (F19), 

( c 2°2) 4 < 1 


n 2(^ s )| 2 | fi 2”(^ s )| 2 D min( a >fl 


(0 £ v £ 1) (F20) 


Analytical bounds for this function are not as easy to determine. For p = 0, 
(v, s ) = |jp2 ( y > s )| * and 


c 2 g 2 


p=0 


= \^~ v 

a* 


or 


(F21) 


where g^, vj has been investigated and tabulated by Case, De Hoffmann, and Placzek 

(ref. 25). They show that gf^-,i'J occurs at v = 0 for a < 1^1% whereas for 

/ |max 

a > 7t2/b it occurs for v between 0 and 1. For a very large, they have 
g max — 4 o;2/ 7 t 2 . The present geometric interpretation (ref. 26) is consistent with all 
these characteristics. The radius of curvature of the curve C 2 given by equation (F17) 
is 7 t2/b at (a',/3') = (0,0). For a' very large, /3 T — v/2; thus, the minimum squared 
distance from (a,0) to (a',p') approaches u-2/4, in agreement with equation (F21) 

3- n< 4 g max ““ 
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Note that the exponential factor equation (F8) requires a > 0 in 

order for B^(a,/3) to be bounded since both v and A are nonnegative. On using 
estimates F(14), (F15), and (F20), one obtains from equation (F7) a bound for b2(o!,/ 3) 
which is denoted as B ma J(o!,/3): 


B 


2 

max 


e -4o?A {] 

“( ct R * 0 - -of + 0 2 ( a R - l) 2 }(“ 2 + 0 2 ) 

4D min 4 («,0) 

(a + a D - cr R )2 + P 2 


(or,/3) =< 


(<x > max(l, -ct d + a R )) 


p-4aA 

cfl + j3 2 ^ 

[[“(°R - X ) - ct d] 2 + 0 2 (°R - l) 2 ]E“ - « 2 + P 2 ] 

4D min 4 (<*,0)| 

(a - l) 2 + 

j P 2 

v. 


(p ^ 0; 0 < a < max(l, -a D 


(F22) 


O 

It can be seen that B ma £ depends not only on a and /3 but also on the nondimen- 
sional material parameters Ojj and as well as the slab thickness parameter A. 

The estimate (F22) for B ma: | is not bounded for the following regions in the s-plane: 


Re(s) < < 0) ^ 

s e C 2 > (F23) 

s e Branch cuts of ^(s) u ^02 ( s )_. 

These regions must be handled separately. Even for the general case, where s does 
not belong to any of the regions (F23), it appears difficult to say whether the eigenvalue 
problem has nontrivial solutions. One suspects that it has only trivial solutions for such 
regions since that is the result which has been found for certain special cases by others. 
Lehner and Wing (refs. 6 and 7) have shown this result for the bare slab, whereas Lehner 
(ref. 8) and Hintz (ref. 9) have obtained this result for the slab surrounded by pure 
absorbers. One shows that this result is also obtained for the special case A — °°, that 
is, a thick slab, as follows. 

Since the slab-thickness parameter A appears only in the exponential term of 
equation (F22), it is seen that B ma J(o!,/3) can be made as small as one likes as A — °° 
if s does not belong to any of the regions given in regions (F2 3). For |B max (a!,/3)| < 1, 
the Neumann series solution of the inhomogeneous integral (eq. (F3)) converges to a 
unique solution. (See ref. 27, for example.) Fredholm's alternative theorem (ref. 27) 
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then guarantees that the corresponding homogeneous equation, namely, equation (FI), has 
only the trivial solution. Thus, for s e S2 e , the eigenvalue problem has only the trivial 
solution as A — «>, When s e 82^, the unique Neumann series solution of equation (F3) 
must satisfy the additional constraint (eq. (F4)). By using the condition C 2 (^02> s ) = 
one obtains 


e -2(s« 2 )a /v Q2 _ j- p(s) eie(s) jA 


where 


p 2 (s) = 


Re(no 2 ) - l ] 2 +Im 2 (^ 02 ) 
§*e (^02) + l ] 2 + Im 2 (^02) 


0(s) = tan 


-1 


Im (^02) 
Re(v 0 2) " 1_ 


- tan 


-1 


Irn(^ 02 ) 
Re (^02) + 1_ 


Now since Re^Q 2 ^= 0, one has 


(F24) 


(F25) 


m a a -° 


(Re(n 02 )^0) (F26) 


Therefore, the Neumann series solution is seen to converge to zero as A — «° when 
Re^C^) r 0. Note that = 0 is the branch cut of nQ2(s) which is one of the 

regions given by region (F23) which must be considered separately. When s e C2, 

Uq 2 = V) 0 = 77 = 1 so that p^(s) of equation (F25) becomes 

P 2 (s) = jj 2 2 1 (F27) 

and p = 1 occurs only at s = -a 2 (that is, (a,/3) = (0,0)). One uses equation (F27) in 
equation (F5) and, on taking the limit A — finds that E2 ± (m,s) — 0 for s e C 2 , 
s ^ -a 2 . 

Summarizing the results then for A — °° indicates that the eigenvalue problem 
has only the trivial solution for Re(s) > -a 2 unless s belongs to either the branch cut 
of t'Qi (s) or (s). In order to determine what happens on these cuts, one must write 
equations (63) and (64) with I m± = 0 in terms of the X m (-z,s) functions rather than in 
terms of the Xq(-z,s) function. This will be done in appendix G. When A is not 
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large, others (refs. 6 to 9) have shown for special cases of the present problem that if 
the eigenvalue problem has nontrivial discrete solutions, such conditions occur on the 
real s-axis. For the bare slab, it was shown in references 12 and 13 that these solutions 
lie on the branch cut of I'q(s). In view of these results, it is assumed that the eigen- 
value problem has nontrivial solutions for Re(s) > -02 only if s belongs to either the 
branch cut of z'oi(s) or ^02( s )* 


83 



APPENDIX G 


SOLUTION OF THE ASSOCIATED EIGENVALUE PROBLEM FOR s e S 2i 

In this appendix, the solutions of equations (63) and (64) with I m± s 0 for s on 
the branch cuts of fg m (s) are examined. It is convenient to use coefficients related 
to the original expansion coefficients A m± (t') and a m± (the bar indicates that one is 
considering the associated eigenvalue problem). Recall that the E m± are related to 
such coefficients by equations (50), (E4), and (E5). It was also noted in appendix F that 
the coefficients can be determined only to within an arbitrary factor independent of v. 
By following references 12 and 13, one introduces coefficients B m± as 


A m± (v) = a 2± B m± (i/) 
H± = &2±*>1± 


(Gl) 


9 

The estimate B max (a,/3) (eq. (F22)) was not bounded on the branch cuts of ^Om^)* In 
that estimate the Xq(-z,s) function was used so that the behavior for s inside, on, and 
outside the curve C 2 could be seen. To investigate what happens on the branch cuts of 
^ m (s), one should use the X m (z,s) functions (appendix A) which do not contain the 
branch cuts of ^ m . Also, when z'gi(s) becomes pure imaginary (that is, on its branch 
cut), one cannot include the contribution (the pole at z' = in the integral over the 
contour C 1 of equation (64). Recall that the material properties c m and a m deter- 
mine where on the real s-axis the branch cuts of ^o m (s) lie. The only restriction 
which has been made is that c^ < 1 and this restriction alone does not specify the over- 
lapping of the cuts. It does, however, guarantee that the branch cut of fQi(s) l ies 
entirely to the left of s = 0. 

Consider s e S 2i n S-j^ first. When the relationships (Gl), (50), (E4), and (E5) 
are used in equations (F3) and (F4), one obtains after some algebra and use of the 
X-identities of appendix A, equations for B 2± (p.) and the additional constraint, namely, 
equations (77) and (81). Recall that equations (F3) and (F4) were obtained from equa- 
tions (63) with I 2± = 0. Equations for B^ ± (-p.) and bj ± are obtained in a similar 
manner from equation (64) with Ij_ ± s 0 when the contour C' is collapsed on to the 
interval (0,1) of the branch cut of £2 2 (z',s). These equations are given as equa- 
tions (78) and (79). The normal-mode expansion of the solution of the associated eigen- 
value problem is given in terms of the B m± coefficients by equation (76). Note that 
equation (81) is the exact eigenvalue condition since all material properties have been 
assumed to be known. It determines the values of s, {%}, for which the eigenvalue 
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problem has nontrivial solutions. When s belongs to the branch cut of Vq^, equa- 
tion (81) takes on different values above and below the vqi cut. Therefore, it is con- 
cluded that the eigenvalue problem has only the trivial solution on the branch cut of i'qj. 
On that part of the branch cut of vq2 which is not also a part of the cut, equa- 

tions (77) to (81) require that the limiting values of the coefficients above and below the 
Vq 2 cut be related as 


[ B m±^)] “ ± [ B m±(/ i )] 

[ E l±] + = *[?!*] 


> (G2) 

( Re M = Im M = °) 


that is, where s is real and is given by max£- 02 , -cr i(l " c i)] < s < - 02(1 - C 2 ). It then 
follows from equations (76) and (G2) that the limiting values of 1J> ± (x,/j.,s) for the same 
region are given by equation (96). From the results of references 12 and 13 for the bare 
slab, it is expected that the eigenvalue problem has nontrivial solutions only at isolated 
points {s n ^> which lie on the branch cut of Pq 2 but not on the branch cut of Vqi- 


In the limit C2<?2 a - — 00 which was discussed in appendix F, one sees that equa- 
tion (77) gives B 2 ± (/i) — 0 whereas equation (81), the eigenvalue condition, becomes 


*2(-to») e ( SW 2) a A>2 

X l( -l '02’ s ) ‘'Ol + v 02 + x l( t '02 J s ) ’'oi " ^02 

^c 2 a 2 a - °° and Re (1^2) = Im(^ 01 ) = 0) (G 3 ) 


Equation (G3) is the "thick-slab" eigenvalue condition and for the region of the s-plane 
where it is valid, it can be seen that one obtains an even eigenvalue s n if 


Im 


X 2 (-02,3) 

X l( I/ 02’ s ) ^01 “ ^02 


= 0 


(G4a) 


and an odd eigenvalue s n if 


Re 


x 2( v 02> 

X l{ v 02> 


f) 

s ) 


js n +<j 2 y/p 02 

V 01 ' v 02 


= 0 


(c2 <J 2 ?i ~ °°) (G4b) 
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Note that equation (G3) has the same form as the zero-order approximation of the critical 
condition given in reference 2, except that here one has both even and odd parity solu- 
tions. Numerical solutions of the eigenvalue conditions are discussed in appendix J. 

In the region of the s-plane, s e S^ e D S 2 i5 one is specifically interested in the 
solution on the branch cut of ^(s) which lies to the left of s = -oq, that is, for s 
real and -02 < s ^ -oq. For such values of s, the solution (76) outside the slab is not 
bounded as x — °°, since 


s) = <pi v (fi, s) 




(0 £ v * 1) (G5) 


In addition, the restriction Re(s) > -O’min on both inhomogeneous terms I m± (see 
eqs. (65) and (66)) also indicates that one cannot deform the inversion contour to the left 
of Re(s) = -oq in general. However, when one is looking for the solution inside the 
slab, |x| ^ a, perhaps the inversion contour can be deformed to the left of Re(s) = -oq 
for special values of material properties and/or initial data. For s e S 2 i D Sj e , expan- 
sion coefficients for the solution inside the slab are obtained as equations (82) and (83). 
Note that equation (83) is exactly equation (81) with Xq^(z,s) replacing (Vqi - z)Xq(z,s). 
Recall from equation (A10) that these are the X-functions which are continuous as s — Cq. 
Under the same replacement of Xgq(z,s) with - z)Xq(z,s), equation (82) reduces 

to the equation from which equation (77) was obtained. Equation (83), which corresponds 
to the eigenvalue condition equation (81), determines the pseudo-eigenvalues, that is, the 
values of s, -03 < s < -oq, where \p2±( x >l 1 ) s ) has nontrivial solutions. 
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FORM OF i// ± (x,/i,s) ON THE BRANCH CUTS OF i/ 0m (s) 

In this appendix, the transformed solution xf/ ± (x,ti,s) is put in a form where one 
can see how it behaves on the branch cuts of I'om^)' ( - >ne ex P ec ts that if/ ± contain the 
branch cut of ^oi( s ) since only one of the two discrete modes appears for |x|>a. 

Such branch cuts appeared in the half-space problems solved in references 14 to 16. One 
also expects that the branch cut of ^(s) does not appear in \j/ ± but instead one should 
find poles at s = s n , the places where the associated eigenvalue problem has nontrivial 
solutions. This condition was found for the bare-slab problem analyzed in references 6, 
7, 12, and 13. 

It is not obvious from the equations which determine the expansion coefficients 
implicitly how one should group the terms to show what is expected. Consider first 
i^2 ± (x,/ii,s). From equations (32) and (44), 

<// 2 ± (x,m, s ) = a. 2 ± ^Pv 02 (x,ii,s) ± ^-i/ 02 (x,M, s)J 

+ J a 2±(^ |^2y( x ^’ s ) ± ^2(-zv)( x ’^ s )] dv 

+ F 2±( x,t, 02» s ) ± F 2±{~ x > v 02> s ) ^^ 02 ( x ^,s) 

+ |F2±( X ^> S ) ± F 2±(" x > y > s ) ^2(-y)( x ^ 5 s )] dl/ ( H1 ) 

Note that this equation can be obtained from equations (70) and (72). It is readily shown 
from the definition of the F m± functions (eqs. (46)) and the properties of the C m± 

(eqs. (Bll)) that 


f 2±( x > z/ 02> s ) = F 2±( a ’ z '02’ s ) ± f 2±(" x >-^ 02> s ) 

> 

F 2±{ x ’- V 02’ s ) = F 2±( a ’ - ^02> s ) =•= F 2±{- X > v 02’ s ) ^ 


(H2) 


It follows from equations (H2) that the two coefficients in equation (HI) can be written as 
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■^2±( x > t '02’ s ) = \ j^ P 2±( x ’ t ’02’ s ) ± F 2 ±(- x ’^ 02’ s )_ + f F 2±( a ^02> s ) 

±F 2±( _x » l/ 02’ s ) = \ | F 2±( X >" I/ 02 ,S ) ± F 2±( -x J I/ 02’ s )j ± | F 2±( a ^02, s ) 


where equation (B12) has been used to replace F 2± (a,-yQ 2 ,s). Equation (HI) then 
becomes 

^ 2± (x,/i,s) =[a 2± + | F 2±( a ^02’ s )]^ 0 2 (x ^’ s) * ’ // -^ 0 2 (X,M,S) ] 

+ J Q A 2±(^ fyzvfaw) ± ^ / 2(-z')( x ’ /i ’ s )] dl ' 

+ \ [ F 2±( x ^ 02’ s ) * F 2 ±( -x,_i; 02’ s ) % 2 (x ’^ s) 


+ | [^2±( x ’“ t ’02’ s ) ± F 2 ±(“ x ^ 02> s )] 

+ J q [f 2± (x,z/,s) ^ 2 „(x,jLi,s) ± F 2 ± (-x,v,s) 1 #' 2 (_„)fc,M,s)]dJ> (H‘ 

When s e Branch cut of yQ 2 (s), then z^q 2 = i | ^02 1 ^ or ^ m ( s ) = 0" an< ^ V Q2 = “^ | ^02 1 
for Im(s) = 0 + . Therefore, on going from below to above the branch cut of z^q 2 , it is 
evident that the third and fourth terms in the right-hand side of equation (H4) simply 
interchange whereas those containing F 2± (x,v,s) and F 2± (-x, zqs) are unaffected 

since these functions do not depend on fg 2 . The coefficient of J^ a 2 ± + ^ F 2 ±( a 5 l '02’ s ) 

however changes sign for odd-parity solutions and the behavior of A 2± (f) is not clear 
yet. By comparing equation (H4) with equation (76) for |x| < a, one suspects that 

\*-2± + \ F 2±{*’ v 02> is the coefficient which excites the associated eigensolution 

^ ± (x,/x,s). This is the information needed to group terms in the implicit equations for 
the expansion coefficients. 

Now examine the equations which determine the expansion coefficients. From 
equation (52) upon using the X-identities, the definition of the h m functions (eqs. (80)), 
and the relationship between the E m± and the original expansion coefficients, the fol- 
lowing equation is obtained: 
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0 = 


a 2± + f F 2±( a >"02> s ) 


h 2^02) h 2(~^02) 


V Q1 + v 02 y 01 " v 02 


+ | F 2±( a ^02> s ) 


^2(^02) _ h 2(~^02) 


^01 + u 02 ^01 “ v 02 


£ ^2±(M) + F 2±( a >M,sj] h 2 (ju.) 


dju. 


pi [1 + Vn-i 

± \ F 1± (-a,fz,s) h i(jx)— j dfA * F l±( _a » v 0l* s ) h l(*^)l)- 

0 - ^Q2 f/ 


2^oi 


01 ‘ v 02 


M + ^01 
(H5) 


By following the same procedure with equation (51), one can obtain after making use of 
equation (H5): 


A 2±^) =^§ F l±(- a ^’ s ) 


± k ^2(°°> s ) u 02 " 1/2 h 2W 


2 %Ks) ^2 


r i„ / N 

y 2 N^M\p± + 2 F 2±{ a -> v 02’ s ) 


^2(^02) h 2(-^02) 


u+v 02 v- Vq 2 


* f F 2± (a,» 02 ,s) T + C |a 2± (m) * F 2± (a,n,s)] ha&i) ^ 


pi ^m 2 " 2 9»9 >»s) 

£ F 1± (-w) h l( ,) _ ^ ^ da 


(H6) 


From equation (54), there is obtained: 


Th l(“^0l) a l± “ a 2± + ^ f 2 ±( a » ^0 2 » s )j j^2 (^0 2 ) ± h 2( _z/ 02)j 

+ |F 2± (a,^ 02 ,s) jh 2 (v 02 ) + h 2 (-^ 0 2)] 


(^02 " v 0l) £ \^2±M + F 2 ±( a 5M> s )j h 2 (ft) 

U jl V, 

{ u 02 - v 0 l) £ F l±(- a >M,s) h x (M) - -^ - 2 


F l±(- a ^0D s ) h l(-^0l) ±h l(^0l)] 


(H7) 
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Finally, from equation (53), 

. . . c 2 ct 2 a , s ( a l- a 2) z / v 

a 1±(-v)*1^ a 2±( v ) ^ II 


= ± 


k b) hjiv) 

2 ^(oo, s) N 1 (y) 


\ < [^2± + | F 2±( a ^02^ 


^2(^02) ± ^2(2^02} 
02 v + v 02 


| F 2±( a > y 02’ s ) 


^2(^02) h 2(~^02) 


v ~ v 02 v + 


v 02 


pi r— “I 2(jP^,,(p.,S) 

+ j |A 2± (fi) + F 2± (a,ji,s)| h 2 (m) cTaTv dfX 


C 1 u i 

\ F l±(-a,M,s) h^pi) - 

J 0 v. 


'of - M d M 

, 2 - u 2 M ^ 
02 ^ 1 


F l±(-a,y,s) - ^^F 2± (a,i/,s) e^ 1 / 


(H8) 


By following a procedure similar to that of reference 13, the expansion coefficients 
A. m± (ji) and a^ ± can be written in the form: 


A m±(l 1 ) " 


a 2± + f F 2±( a ^02> s ) A m±W + B m±^) 




a l± = 


a 2± + | F 2±( a ^02^ s ) a l± + bi±(M) 




(H9) 


When equations (H9) are used in equations (H6) to (H8), it follows that 






m± K 


a i± = b i± 


(H10) 


where B m± (p) and b^ ± are the expansion coefficients of the associated eigenvalue 
problem given by equations (77) to (79). The coefficients B m± (|U.) and bi ± are found 
to be given by equations (85) to (87). The coefficient a 2± is obtained from equation (H5) 
ds 
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a 2i - f Ma,, 02 ,sj| 1 + f B 2± M h 


dp. 


2^ [i + v Q1 


= -| F 2±( a ' 1 '02> s ) %i ( ?4 T ^Ol'- ° ol] ‘ CM + F 2±( a ."> s )) h 2(M) ^ 

pl JUt "H Z^n 1 2 Z^n 

T j 0 F 1± (-a,p,s) h^p) — — ^ dp + F 1± (-a,i/ 01 ,s) h 1 (i/ Q1 ) - 


„ 2 _ v 2 
"01 02 


01 

(Hll) 


It can be seen from equation (81) that the coefficient of 


a 2± + | F 2±( a ^02 5 s ) 


in equa- 


tion (Hll) is the eigenvalue condition, and it will be zero at the places where the associ- 
ated eigenvalue problem has nontrivial solutions. Equation (Hll) is the same as equa- 
tion (88). The solutions \p m± (x,ii,s) can now be written as equations (91) and (92). 
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BEHAVIOR OF ip ± (x,p,s) ON INVERSION CONTOURS 


In this appendix, several points concerning the behavior of i^ ± (x,p,s) on the inte- 
gration contour of the inverse Laplace transformation and some parts of related deformed 
contours are discussed. First, one looks at the behavior of ^ ± (x,p.,s) as |s| — °° with 
Re(s) = y, a large finite positive number. It will be seen that ^ ± (x,/r,s) is not neces- 
sarily 0^|). Such parts of i// ± (x,p,s) are inverted separately and the resulting solutions 
are shown to satisfy the uncollided transport equation. Then one considers how !// ± (x,/a,s) 
minus the uncollided term ^ u± (x,p, s) can be deformed around the poles and the branch 
cut of ip ± . 

It is of interest to examine the behavior of ify ± on the contour Re(s) = y as 
| s| — where y is finite (see fig. 4). For such cases, s e S^ e n S 2 e and the solu- 
tions i// m± can be seen from equations (32), (39), (44), and (45) to be 


^ 2± (x,fi,s) 


Jq |^2±W + F 2±( x ’ y ’ s 3 e ^ 
± < f 0 ^2±( I/ ) + F 2± (-x,^s)] ^ I 


<P 2 j,(m, s) dv 


V 

VZvi-V’S) dv 


(ID 


and, for x > a, 


\p 1± (x,i±,s) =±y jA 1± (-i/) - F ± (-a,^,s) + F 1± (-x,-^,s)j e ^ $ I <P\ v {p,a) d » 

pl (s+o^x/v 

± j F 1± (-x,v,s) e v '' <p 1 v (-il,b) dv 


(12) 


with an equation similar to equation (12) for x < -a. One sees then that the coeffi- 
cients A^ ± (-v), A 2 ± (z'), and the F m± functions are needed. The expansion coeffi- 
cients are given implicitly in terms of the F m± as 
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/il 2± (fi,s) -(s+c 2 )a jv. 




Oj(«,s) uX 0 (-m,s» e ‘( St< ’2)7 1 ' f l -(s*a 2 )a/i’ 

S .M IS i ivtw. d - 


c 2 a 2 fil(°°, s ) 


N 2 (m) 


( 13 ) 


and 


a i±(-m) 


(s+CTj) a Ai 

fj.I 1± (fx,s) e 7 cxax N 2 Qx) "( CT 2 _a l) a / (X 


A 2 ± (m) 


Nx(m) * c 2 ct 2 Nx(m) 

(s+o^a/n ! A M x e -( S+<T 2) a / 1 ' 

Xo(-n.s) N^y J 0 5^2 


fxk e 


(14) 


where the I m± (/i,s) are given by 


Coffn (s+CT X )a/| 

W^®) = c^r N i ( ^ F i±(- a ^> s ) eV ' 


v- 


^§Ki| X 0 ( ' M ’ S) Iq F 2± (a>I/ ’ S) 6 Xo(-^s) <p 2v {-n, S) di/ 

r 1 F l±(“ a > I; > s ) e( S+CTl ) a A 
x o (_ ^’ s) J 0 x 0 ( ~ ;s ) *i^ 8) ^ 


c l a l 


(15) 


and 


Il±(M,s) F 1± (-a,/a,s) e F 2± (a,M,s) e ^ /* 


V- 


+ iylo F2±(£W) e (s+<72)ay/ " x o ( -^ s) ^ 2 >»®) 

(s+a^a/12 

k %(« s) x pi F 1± (-a,v,s) e' 7 ' 


dv 


c l u l ^2(“ s ) x 0( 


i_f l 

-ju,s) J 0 


X 0 (-y,s) 


<Pjj,(-M-,s) dv 


(16) 
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The behavior of various functions which appear in equations (II) to (16) as |s| — 
Re(s) = y is 


O m (z,s), A m (^,s), k - O(s) 

X 0 (z,s) - 0(1) 
and for 0 S n,v = 1, 

<P mv (lJ.,s) - 0(s) ^ 

> (17) 

<P mZ /(-M,s) - 0(1) 

The F m± functions appear with an exponential factor and the combined behavior in the 
same limit is 


F m±( x >^ s ) 




On using equations (17) and (18) in equations (15) and (16), one finds that 


Ml2±(M,s) 

n 2 (m) 


F 1± (-a,p,s) 


e (s+0])a/M 



(19) 


and 


m.i 1 ± (m,s) 

N x (p) 


+Fi ± (- a ,ja,s) 


e -( s+<, i) a A 


± 


F 2± (a,M,s) 




(HO) 
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The coefficient A 2 ± (/i) is obtained from the integral equation (13). Since the kernel of 
this equation is also the first term of the Neumann series solution will give the 

behavior of A 2 ± (ju) as |s | — It follows then from equations (13) and (19) that 


A 2 ± (m) 


F 1± (-a,M,s) 




-o(i 


( 111 ) 


Using equations (110) and (111) in equation (14) yields 

r ~i “( cr 2“ cr l) a /^ -2(a 2 -a 1 )a//i 

[A 1± (-m) ± F 1± (-a,jtx,s)] - ±F 2± (a,M,s) e ' ± F 1± (-a,/i,s) e ' 


- off 


( 112 ) 


Equations (111) to (112) are used in equations (II) and (12) to get 


r 


i// 2± (x,/x,s) - < 


F 2 ± (x,m,s) + F 1± (-a,M,s) e 


-( a 2- a l) a /M 


-(s+Cr 2 )x/|U. 

(s + cr 2 ) e ' 


~(s+a 2 )x/ii ( cr 2“°'l) a / ( 

+ cr 2 ) e ' F 2 ± (-x,-ii,s) + F 1± (-a,-M,s) e x 7 / 


±(s 

v 




(M > 0) 

(M < 0) 

(113) 


and 


(s + ^l) 




F 1 ± (-x,-m,s) - F 1± (-a,-M,s) 


^i±(x,M,s) - / 


* F 2i (a,, )S ) * Fj ±( -a,fi,s) 


b(s + 


'( s+a 1 ) x y/ i 


F 1± (-x,-/i,s) 


(jLt > 0) 

(M < 0) 

(114) 


when x > a. For x < -a, has a similar form. Upon using equation (18) for the 

F m± functions, one finds that equations (113) and (114) can be written as equations (98) to 
(101) where the symmetry properties of f m± (x,/z) have been used. It can be seen from 
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equations (98) to (101) that if/ m± (x,ii,s) is not necessarily In fact, if f m± (xo,/i) 

contains 5{x - xq), then i// ± (x,/i,s) is 0(1) as |s| — °°, Re(s) = y. The parts of 
4/ ± (x,ii,s) which are not can be inverted by inspection after a change of variables 

is made. 

For x > a, p. > 0, and all s, <// u± (x,p,s) is defined as 


nx -(s+ai)(x-xo)Xi , 

!// u± (x,m,s) ® - j e / fi ± (x 0 ,M)/dx 0 

ci 

" ( a 2 " °l) ( a ” X )A pa -(s+a 2 )(x-x 0 )/p. 

+ p J e ^ ' f 2±( x 0^) ^0 

e -(a 2 -a 1 )2a/M -(s+<t 1 )(x-x 0 )/m 

+ J e 11 fi ± (x 0 ,M)dx 0 


(115) 


which gives the parts of equation (100) which are not for |s| — °°, Re(s) = y. Now 

the change of variables 

x - xq = lit (116) 

where t = 0 since x = Xq and p > 0 is made. Equation (115) then becomes 

r (x-a)/p -(s+aV)t 

^ u± (x,/r,s) = J e v 1 f 1± (x - lit, jll) dt 

r (x+a)//i -(s+ojjt _ ( a 2 _ff 1 ) (a _ x)/u 
+ \ e e f 2± (x- lit, ii) dt 

J (x-a)//i 

-fs+a-t'it -(ao-CTi^a/u 

+ \ e V ' e ^ i; ' f 1± (x - /it, ju.) dt (117) 

J (x+a)/u. 

which is easily seen to be equation (103) with ^ u± (x,q.,t) given by equation (106). For 

M < 0, 

Xq - X = |/i|t (118) 

is used. It is seen then that the results given as equations (103) to (107) follow. 


96 



APPENDIX I - Continued 


Another point to be discussed in this appendix is the contribution to equation (114) 
from the contour Cp (see fig. 4) around the right-hand end of the branch cut of ^oi(s) 
as the radius p goes to zero. This branch point is located at s = -ct-j^I - c^; thus, 

s + a-^l - Cj) = pe icp (119) 


Here ^(s) — °° as p — 0 as 


2 __ c l g l 1 
01 o 3 pe i (p 

The branch cut has already been picked so that z'gj(s) is real when 
greater than -cr-j^l - c-^). The integral 

2 ^ J ^ ± (x,p,s) e st ds = p |i/ ± (x,M,s) e s ^j e 1(p dcp 

Cp -7T 

with s given by equation (119), is zero in the limit p — 0 if 


(-TT < cp < n) (120) 
s is real and 


( 121 ) 


lirn [p^ ± (x,p,s)J =0 


(122) 


independent of cp. As pointed out in the text, the point s = -Oj(l - c-j^ may happen to 
satisfy the eigenvalue condition (eq. (81)). One assumes for the moment that it does not 
and shows later what changes are required if it does. The function f2j(<»,s) — 0 as 
p - 0 as 


%(°°,s) pe W 


(123) 


so that 


v 0l n l(°°> s ) 


S = -CT 1 (l-Cl) 


c l g l 

3 


(124) 
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At this branch point s e S-q n S 2 i so one needs to show the behavior of many functions 
given in the text as p — 0. This behavior can be given in terms of the behavior of 
and In the relationships which follow, quantities which are functions of s will 

be given as O^qj^, O(l/^0l)’ 0(%)» 0(1), etc., as s cr^l - cj). For example, 

f2^(°°,s) Vq^ — Finite — 0(1) (27) 


where the given equation number is that from which the relationship can be seen. 


X m(-M,s) - 0(1) 


x m( ±zy 02» s ) - 0(D > 

x m( ±l/ 01’ s ) - °( 1 Aoi) > 


(A8a) 


k - 0(1) 


(C7) 


h 2 (w) - 0(1) 

N m (pt) - 0(1) 
hi 00 - o(«!) 


(w = ±v Q2 and m(O^ju^I)) 

> (80) 


h i(±^oi) * °(iAoi) 


b 2± (i H) - 0(1) 


(77) 


b i± - o^oi) 


(79) 
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B 1± (-m) - 0(1) 

1^2 ± (x,M,s) - 0(1) 

> 

^ 1± (x,/i,s) - O(v 01 ) 

J 

F m±( x >M,s) - 0(1) 

F 2±( x ’ ±l/ 02’ s ) ®(1) 
F 1± (x,±y 0 i,s) - O(l/y 0 i) 

B 2± {v) - 0(1) 


B 1± (-y) - 0(1) 
a 1± - 0(1) 

> 

a 2 ± ~ 0 ( 1 ) 

-/ 

h± - °d) 

> 

h± - o(D 

J 

F ±( -a ’^01’ s ) - O(lAoi) 

F ± (-a,y,s) - 0(1) 


(78) 

(76) 

> (46) 

(85) 

( 86 ) 

(89) 

(90) 
(46) 
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[ b l± " F ± (-a^ 0 l,s)| - O(V 01 ) 


a 2± + | F 2±( a ’ y 02 5 s ) 


- 0 ( 1 ) 


(87) 


( 88 ) 


i// 2± (x,m,s) - 0(1) (91) 

\js 1± {x,i±,s) - O(^ 01 ) (92) 


From these last two relationships and equation (124), it follows that 


pi// ± (x,p,s) - {p 0(1) 


(125) 


so that equation (122) is satisfied. Therefore, there is no contribution from the inte- 
gral (121) for the case when s = -0^(1 - cj) does not satisfy the eigenvalue condition (81). 

If the point s = -cr-^1 - CjJ happens to satisfy the eigenvalue condition, the denomi- 
nator of ja 2 ± + |’ F 2±( a ’ y 02’ s )j is equivalent to the eigenvalue condition (81) van- 

ishes. It can be seen from equation (88) that the limiting form of this condition at the 
branch point is ai ± = 0 and additional discussion about this condition is given in 
appendixes J and K. If one considers, for such cases, the function 


i// ± (x,p,s) 


a 2± + | F 2±( a ^02^ s ) 


i//,(x,p, s) 


(126) 


instead of i// ± (x,ju,s) as the integrand of the integral (121), it follows that in the limit 
p — 0, the contribution from such an integral vanishes. The part which has been sub- 
tracted from ^ ± (x,p,s) in function (126) is considered separately and would appear to 


have a pole because of the zero in the denominator of 



bution therefore does not vanish in the limit p — 0; in fact, its contribution looks like a 
discrete residue term. However, the point is not isolated ^remember that this is the 
branch point of Vqi at s = -cr^l - so it should be understood that its contribution 
is included in the branch-cut integral term of equation (124). The numerical results 
indicate that such points occur. 
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APPENDIX J 


EQUATIONS AND PROCEDURES FOR COMPUTATION OF EIGENVALUES 


The equations from which the time eigenvalues {s n } are determined for 
s e S;q fi S 2 i are equations (77) and (81). When s e Sj e n S 2 i, the corresponding equa- 
tions are equations (82) and (83) and they determine what has been called the pseudo- 
eigenvalues. These equations are solved numerically by using the procedure of refer- 
ences 12 and 13. As stated previously, these equations can be written in terms of the 
nondimensional quantities introduced in equations (118). By making the substitution 


b ± (m) 


B 2± (/a) h 2 (p) fa i(^Q2, s ) -(f) 

( y 01 + ^02 f fi l(°°> s ) 


(Jl) 


it follows that equation (77) can be written for £ real and max(-crj) + Or, 0) < £ < 1 
(that is, on that part of the branch cut of vq2 which is not also part of the branch cut 
of ^qi) as 

B ± (p) = -g(ju) 


^ g± * 1 ^02l g T 


p 2 + 


"02 1 


If 1 

2 Jn 


B+(^) 


dv 


v + \i 


(J2) 


where 


g (m) 


g± = 


%( I/ 02^) 

x 2 (-m,0" 

%(°°,C) 

Xi(-M,C)_ 


jig - Q e~ 2 ^ M 2 + |^q 2 1 2 
n 2~(^ (M + ^Ql) 2 


l^l(i/ 02 ,C) Im 

x: 2( I/ 02^] 

e ^ A A02 

| «l(°°,C) Re 

Xl(^02^J 

1 ^01 “ v 02 






Now one defines 



as 



d.v 


which is the eigenvalue condition (eq. (81)) if 



= 0. 


(J3) 


(J4) 
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Equations (J2) are reduced to two sets (±) of N equations in the N unknowns 
and B_(ju.i) where i = 1,. . .,N (see, for example, ref. 28) given by 


B ±( M 0 = 




+ho2|‘ 


if K.^1 

2 l-j i M-; + IJ'i 
i=l j 


(J5) 


where Rj are the weighting functions for the numerical integration scheme which is 
used. Equation (J4) is written as 

N 


3=1 


(J6) 


Since a search must be made for values of £ for which A^ ± = 0, one must be able to 
compute all quantities which appear in equations (J3) for any value of £ in the range 
given by equation (120)). These quantities are computed as follows: 


^l(^02>g) _ g + g D ~ £°R 
^l(°°iO C + " a R 


^ 2(^,0 


c 2 a 2 



(J8) 


The functions yQ m are determined by ^ m (^0m^) = 0 and they are computed numeri- 
cally by using the Newton -Raph son iteration (ref. 28) on the nonlinear equations 




(J9) 


and 


vqI tanh 


-1 


1 

V 01 


C + a D 
°R 


(J10) 


The X-functions are computed from equation (A8a); namely, 
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For £ 


where 


For £ 


where 


x m< z >«=rb exp ^ lo log ‘ 



d^ 

«m(^C)_ 

1 / - z 


(Jll) 


and z real where z = -ji, 0 = /i = 1, one has from equation (Jll) that 




(J12) 


e 2 (>Z,v) = tan 1 ( / n 

\£ - v tanh 1 yy 


0 l(^ a R^ a D^) = tan " 1 ( 


R 


irv/2 


£ + °D “ a R zp tanh 1 vj 


(J13) 


real and z = v Q 2 , one calculates the real and imaginary parts of 


Re 


Im 


x 2( l/ 02^) 

_ e r l/» 

"\ 

r 2 
PAC Z. 

X ir02’0_ 


tUQ 

77 

> 

X 2^02’0 

- ^ 

. r 2 
qi n _ 

x l(P02.C). 


oJLil 

71 

J 


X 2(PQ2 ■?) 

x l( 1 '02.?) 


as 


(J14) 


r !-£ 


y 02| 2 


r 2 = |^02 1 “ 0 l(^°R’ a D^) ^ 


dv 




\ v 02\‘ 


(J15) 
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Integrals in equations (J12) and (J15) are computed as 


1 M 

J dv =J J R ife(^i) ' 0 l(^ a R’ a D^i)l^i) (J 16 ) 

° L i=l L J 


where Rj are again the weighting functions for the numerical integration scheme. 

In all numerical integrations, Gauss’ method (ref. 28) was used. For integrations 
in equations (J5) and (J6), the interval (0,1) was split into four intervals 


(0,1) = (0,0.05) + (0.05,0.1) + (0.1, 0.9) + (0.9, 1.0) 


(J17) 


and a 10 -point Gaussian formula was used in each subinterval. For integrations in 
equations (J12) and (J15), the interval (0,1) was divided as 


(0,1) = (0, 0.1) + (0.1, 0.9) + (0.9, 0.99) + (0.99, 0.999) + (0.999, 1.0) 


(J18) 


and in each of these subintervals a 10 -point Gaussian formula was also used. The sub- 
division (J18) is the same as that used by Kowalska (ref. 29) and the X-functions calcu- 
lated here agree with those she gives to all figures which she quotes except for the real 
and imaginary parts of X m ^Q 2 ,£). She apparently used T 2 instead of ^/ti in equa- 
tions (J14) to obtain the numerical values for the real and imaginary parts given in part II 
of reference 29. Since her later published critical-slab results (ref. 24) agree with those 
of Mitsis (ref. 22) for a bare slab, it is expected that this oversight was corrected. 


Conditions (82) and (83) which determine the pseudo-eigenvalues for s e Sj e n S£i 
lead to very similar equations which will not be written down. In this region, the real 
s-axis corresponds to 0 = £ = -ctq and such equations need be considered only if 
-CTj} > 0. 


The procedure used to calculate the eigenvalues £ n is as follows. For fixed 
values of A, cr R , and a D , one selects a number of £ values in the interval given by 
equation (120). For each of these values, one obtains | ^02| anc * ^01 f rom equa- 
tions (J9) and (J10) by iteration (Newton-Raphson). Equations (J13) are evaluated at each 

X 2 (-Mj,g) 


of the 50 Gaussian integration points iq (0 < zq < 1). Next, the 


are 
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APPENDIX J - Continued 


calculated for each of the 40 Gaussian integration points, j u 


(0 < jUj < l) by using 


equation (J16) in equation (J19). The real and imaginary parts of 




are com- 


puted in the same way from equations (J14) to (J16). 


equation (J3) at each of the 40 points 
N equations in the N 


Now g(Mj) can be computed from 




and all the coefficients in the two sets (±) 
(eqs. (J5)) can be evaluated. 


of 


unknowns and B Ax A 

These two sets of simultaneous equations are solved numerically for which are 

then used to compute A^ ± from equation (J6) at the selected values of In this way, 
one locates the zeros of A^ ± approximately. A new set of £ values, located about 
each approximate £ n , is selected and the process is repeated. For the present compu- 
tations, the £ n were located to three figures. Discussion of the computed results is 
given in the text. The calculations were done on a Control Data 6600 computer system 
at the Langley Research Center. 


In appendix G, the thick-slab eigenvalue condition was given as equation (G4). Note 
that g ± quantities given by equations (J3) are, within a factor, exactly the quantities 
needed in equation (G4). Therefore, the thick-slab approximation eigenvalues are 
obtained from 


g± = 0 (J19) 

as would be expected from equation (J4). 

The bare-slab eigenvalues are obtained when or = 0 and it is easily shown that 
in this case, equations (J5) and (J6) no longer depend on o^; that is, for Or = 0, these 
equations do not contain Or. 

It was noted in the text and in appendix I that the branch point of vqi located at 
s = -oj(l - cjJ may happen to satisfy the eigenvalue condition which can be seen from 
equation (88) to be 

a 1± = 0 (J20) 

when Vqi — °°. This point corresponds to £ = -o^ + or and it can be shown that equa- 
tion (J20) then determines values of £ = £ n which depend on neither Oq nor Or; that 
is, if one uses £ = -Oq + Or to eliminate Oq from the condition (J20), Or drops 
out of the equations. Equation (J20) determines the values of £ at which eigenvalues 
disappear into the right end of the branch cut of Vqi- Also note that the limiting form of 
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APPENDIX J - Concluded 


the pseudo-eigenvalue condition for s = -o^, which corresponds to £ = -op, determines 
the values of £ where the pseudo -eigenvalues disappear into the left end of the branch 
cut of z^Ql- Such points, as well as those given by equation (J20), are labeled with an 
asterisk in figures 7 and 9 to 11. 
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APPENDIX K 


REMARKS ON EIGENVALUE— BRANCH-POINT COINCIDENCE 


In this appendix, a few remarks concerning the situation when the eigenvalues (or 
pseudo-eigenvalues) disappear into the branch cut of z'qj are made. This situation is 
somewhat analogous to that encountered by Hintz (ref. 9) for the slab surrounded by pure 
absorbers. He could not say whether a bare-slab eigenvalue (which does not depend on 
Ojq^ that happened to coincide with -Oq belonged to the point spectrum or the continuous 
spectrum for his problem. In the present problem, the eigenvalues coincide with a 
branch point as they disappear into the branch cut of z^. A numerical study of the 
branch-cut integral in equation (124) has not been made nor has the eigenvalue condition 
on another Riemann sheet been investigated. It is suspected that there is no drastic 
change in the shape of the solution given by equation (124) when an eigenvalue disappears 
into the branch cut of Vqi and such studies would resolve this point. It was pointed out 
in appendix J that the condition (J20), which determines whether the point s = - 04(1 - c^, 


(? = _cr D + °r) i- s a zero of the denominator of a. 2 ± + |'F2±( a , i; 02’ s ) 


given by equa- 


tion (88), depends neither on a D nor a R explicitly. In appendix I, it was indicated 
that the contribution from such points should be included in the branch -cut integral since 
it arises from the integration around the branch point. One understands then that such a 
contribution is included in equation (124) if s = -a^(l - cj_) happens to satisfy equa- 
tion (J20). How such zeros of equation (J20) behave or appear in the solution after 
passing through the branch point as the material properties are varied has not been 
studied here. 


If one considered the problem of a finite slab with symmetric reflectors of finite 
thickness, then what is happening at the places where the eigenvalues coincide with 
vqI = °° might be deduced. In such a problem, the solution probably does not contain 
the branch cut of Vq±, but instead has discrete eigenvalues along it. Even though there 
is another parameter in the problem, the reflector thickness, one might be able to carry 
out a numerical study of all the eigenvalues. 
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